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A digital  computer  program  has  been  developed  for  simulating 
the  launch  and  recovery  of  VSTOL  aircraft  operating  from  small 
ships.  The  program,  known  as  VOLAR  (Vought  Launch  and  Recovery 
Dynamics  Program),  employs  the  computational  technique  of  non- 
linear covariance  propagation.  This  permits  the  time  histories 
of  the  means  and  variances  of  all  system  state  variables  to  be 
computed  from  a single  run,  as  opposed  to  the  Monte  Carlo  technique, 
which  requires  multiple  runs,  plus  subsequent  averaging.  Typically, 
VOLAR  requires  approximately  7 percent  of  the  computer  time  required 
for  comparable  Monte  Carlo  simulations. 


The  program  includes  a general  airframe  mathematical  model 
suitable  for  helicopters  or  fixed-wing  aircraft.  The  airwake  of 
a small  ship  is  also  modeled  and  ship  motion  models  are  included. 
Two  alternative  models  of  the  human  pilot  are  supplied,  one  is 
based  on  verbal  adjustment  rules,  the  other  is  based  on  optimal 
control  theory,  utilizing  performance  index  parameters  deduced 
from  manned  simulator  experiments.  The  program  is  demonstrated 
for  AV-8A  recovery  on  a small  ship.  The  trends  predicted  by 
VOLAR  are  shown  to  agree  with  flight  test  data. 


Volume  I of  this  report  describes  the  program,  Volume  II  is 
the  user's  manual . 
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The  reader  who  wishes  to  get  a quick  overview  of  this  report 
without  going  into  details  should  read  Sections  14,  1,  13,  6,  8, 
and  10  in  that  order. 
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1.1  Objective  of  this  Study 

The  study  reported  here  supports  a broad  group  of  related  research  and 
development  efforts  which  are  ultimately  directed  at  enlarging  the  operational 
limits  of  sea-based  VSTOL  aircraft.  Particular  emphasis  is  placed  on  opera- 
tions from  small  ships.  It  is  desired  to  assess  the  feasibility  of  such 
operations  in  quantitative  terms.  For  some  applications  these  terms  describe 
the  probability  of  safe  operation  of  a given  airplane  from  a particular  type 
of  ship  in  a specified  wind  and  sea  state.  Alternatively,  quantitative 
stability  and  control  criteria  may  be  sought  which  yield  a specified 
probability  of  safe  operation.  In  either  case  a mathematical  model  of  the 
dynamics  of  the  launch  and  recovery  system  is  required.  The  principal 
objectives  of  this  report  are  to  develop  such  a system  model  and  to  demonstrate 
its  use. 

1.2  System  Representation 

Figure  1 shows  the  major  dynamic  elements  of  our  mathematical  model  of 
the  VSTOL  launch  and  recovery  system.  The  positional  coordinates  and  velocity 
components  of  interest  are  denoted  as  state  variables.  For  example,  the 
position  of  the  airplane  c.g.  in  inertial  space  at  any  given  time  may  be 
denoted  by  three  state  variables  x,  y,  z,  which  are  conveniently  expressed  as 
elements  of  a column  state  vector  {x,  y,  z}T*.  Other  system  parameters  are 
adjoined  as  necessary.  For  example,  the  ship  inertial  position  and  velocity 
components  xc,  yc,  zc,  x , y , t can  be  adjoined  to  augment  the  state  vector 
to  {x,  y,  z,  x , y , z , X , y , i } . The  number  of  elements  of  the  state 

J d J ) J d 

vector  that  is  required  to  characterize  the  complete  system  is  called  the 
dimensionality  of  the  system.  The  overall  system  model  is  written  in  terms 
of  state  variables  for  convenience  in  programming. 

The  most  significant  inputs  come  from  the  atmospheric  and  sea  motions, 
which  at  any  given  instant  produce  certain  changes  or  perturbations  of  ship 
and  airplane  state  variables.  Other  inputs  are  manifested  as  "noise"  within 

♦Throughout  this  report  column  vectors,  e.g.,  |yj  are  writter  as 
transposed  row  vectors,  e.g.,  {x,  y,  z)T  to  save  space. 
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the  system,  e.g.,  noisy  sensors  of  airplane  state  variables  such  as  angle-of- 
attack  and  sideslip  angle.  For  simplicity,  these  noise  sources  are  not 
explicitly  shown  in  Figure  1.  The  principal  inputs  (wind  and  waves)  are 
correlated  because  the  sea  state  depends  upon  the  steady  wind  and  the  r.m.s. 
magnitude  of  the  wind  fluctuations  (atmospheric  turbulence)  is  approximately 
linearly  proportional  to  the  mean  wind  speed,  as  shown  in  Reference  1.  A 
further  correlation  results  from  the  wind/ship  interaction  which  produces  an 
airwake  which  affects  the  airplane  in  the  vicinity  of  the  ship,  particularly 
when  approaching  from  the  lee  side. 

Mathematical  models  of  each  system  element  are  given  in  subsequent 
sections  of  this  report.  The  block  diagram  of  Figure  1 is  oeneral,  and 
includes  some  elements  which  may  be  omitted  from  some  launch  and  recovery 
systems,  e.g.,  the  LSO  (Landing  Signals  Officer)  and  the  MOVLAS  (Manually 
Operated  Visual  Landing  and  System)  which  permits  the  LSO  to  communicate 
visually  with  the  pilot  in  addition  to  his  verbal  radio  instructions.  Such 
pilot  aids  enlarge  the  portion  of  the  total  system  state  vector  that  can  be 
perceived  by  the  pilot.  This  portion  is  called  the  display  vector.  Its  size 
quantifies  the  complexity  and  sophistication  of  the  display  and  indicates  the 
capability  of  the  pilot  to  sense  his  environment.  Thus,  for  night  approaches 
the  display  vector  is  generally  smaller  than  for  approaches  in  conditions  of 
good  visibility. 

The  mathematical  description  of  the  system  requires  two  further  vectors. 

These  are  firstly  the  control  vector  (u > = (6-j , 62>  6n^’  where  6i’  62* 

6 are  the  aircraft  control  deflections,  and  secondly  the  noise  vector 

" T 

(w>  = {w-j , w2 wm)  where  w^,  w2,  wm,  are  independent  white  noise 

sources.  Physically  all  noise  sources  are  colored,  e.g.,  atmospheric  turbulence, 
sensor  noise,  pilot  "remnant"  (Reference  2).  It  is  customary  to  model 
these  noise  sources  as  outputs  of  hypothetical  filters  forced  by  {w},  instead 
of  the  actual  noise  w'  , where  the  prime  denotes  that  the  noise  is  colored. 

1.  Houbolt,  J.  C.,  Atmospheric  Turbulence,  AIAA  Paper  72-219,  Jan.  1972. 

2.  Kleinman,  D.  L.,  S.  Baron,  and  W.  H.  Levison,  An  Optimal  Control  Model 
of  Human  Response,  Part  I.  Automatica,  Vol.  6,  1970,  p 357-369. 


It  is  frequently  convenient  to  split  the  control  vector  (u)  into  two 
components,  one  derived  from  feedbacks  of  components  of  {x},  the  other 
representing  open- loop  or  pre-programmed  components,  which  can  be  included 
with  {w}.  This  allows  the  system  linearized  dynamical  equations  to  be  written 


{x}  = [F(t)Hx)  + [r(t)]{w} 


(1.1) 


where  (x)  has  been  augmented  to  include  (u>  , and  F includes  additional 
equations,  not  stated  here,  relating  (x)  to  { u)  and  {w1}  to  {w}. 

1.3  Statistical  Descriptors 

Since  some  of  the  inputs  to  the  system  are  random,  meaningful  quantitative 
results  can  only  be  obtained  through  the  use  of  statistical  quantities.  A 
scalar  random  quantity  x is  typically  described  by  its  mean  or  expected  value 
denoted  as  y or  E(x)  and  its  variance,  E (x  - 7)^  . Similarly  a n-component 
random  vector  (x)  may  be  characterized  by  a mean  vector  E{x},  and  a n x n 
matrix  known  as  the  covariance  matrix  P({x}),  defined  as: 


p({*})  = E |(x1  - (xTxl)(Vxn) 


[^n"xn)  (*!-*!>’  * 


(Y*n>2 


(1.2) 


The  principal  diagonal  of  P ( {x >)  defines  the  variances  of  all  the  state 
variables . 

The  mean  and  covariance  can  be  obtained  in  two  distinctly  different  ways, 
as  follows. 

1.4  Calculation  Procedures 

Monte  Carlo  Simulation:  Repeated  runs  are  performed  with  sample 
statistical  inputs,  and  the  results  are  averaged.  Several  hundred  runs  are 
typically  required  to  ensure  adequate  confidence  in  the  averages  used  to  form 
the  means  and  covariance  matrices  as  functions  of  time.  Monte  Carlo  simulation 
imposes  no  restriction  on  system  linearity. 
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Covariance  Propagation:  The  advantage  of  this  method  is  that  it  computes 
the  means  and  the  covariance  matrix  as  functions  of  time  without  the  necessity 
for  running  multiple  cases  and  subsequent  averaging.  Only  two  computations 
need  be  performed,  one  for  the  mean  and  one  for  the  variance.  Thus 
considerable  savings  of  computer  time  are  possible.  Starting  with  an  Initial 
estimate  of  the  covariance  matrix  P at  the  time  tQ,  P at  successive  instants 
is  computed  from  the  "covariance  propagation  equation".  For  simplicity,  this 
is  quoted  below  for  the  case  when  u = 0,  assuming  the  system  is  linear  and  the 
noise  is  comprised  of  a deterministic  component  b,  plus  white  noise  w. 

For  this  case  equation  1.1  can  be  written  as  follows, 

x = Fx  + rw  + rb  (1.3) 

The  parentheses  are  dropped  for  brevity.  F,  r,  are  matrices  and  both  F 

and  r may  be  time-varying  ; x,  w,  b are  vectors  as  previously.  The  covariance 

propagation  equation  is 

P = FP  + PFT  + rQrT  (1.4) 

where  Q is  the  spectral  density  matrix  of  w,  i.e.  E [w(t)  w(t)t]  = Q(t).6(t-r) 
where  6 is  the  Dirac  delta  function. 

By  integrating  equation  1.4  the  time  history  of  P can  be  obtained.  The 
time  history  of  the  mean,  x,  can  be  obtained  by  a modified  form  of  equation  1.3. 

x = F x + r b (1.5) 

Thus  the  desired  statistical  information  [x  (t)  and  P (t)]  can  be  obtained 
by  integrating  Equations  (1.4)  and  (1.5). 

Published  studies  of  the  covariance  propagation  method  have  indicated 
that  it  requires  only  approximately  1/10  to  1/15  of  the  computer  time  required 
for  Monte  Carlo  simulations,  (e.g..  Reference  3).  Covariance  propagation  also 
simplifies  data  handling  and  storage  because  it  works  directly  with  the 
statistical  quantities  that  are  desired,  avoiding  the  computation  of  these 

3.  Gelb,  A.,  and  R.  S.  Warren,  Direct  Statistical  Analysis  of  Nonlinear  Systems: 
CADET,  AIAA  J.,  Vol . 11,  No.  5,  May  1973,  p 689-694. 
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quantities  as  supplementary  items  after  the  numerous  sample  time  histories 
have  been  run. 

Quasilinearization:  Because  Equations  1.4  and  1.5  are  linear  some 
special  adjustments  must  be  made  in  order  to  apply  covariance  propagation  to 
nonlinear  systems  such  as  the  VSTOL  launch  and  recovery  system.  Significant 
ncnlinearities  are  built  in  to  this  system  through  the  thrust  limitations  of 
the  aircraft.  Aerodynamic  nonlinearities  may  also  be  important.  To  explain 
how  nonlinearities  can  be  handled  within  the  framework  of  Equations  1.4  and  1.5 
let  us  consider  a specific  example. 

Figure  2 illustrates  a nonlinear  CL  versus  a variation.  This 
characteristic  is  approximated  by  a purely  linear  gain  which  is  selected  to 
give  the  best  least-squares  approximation  to  the  output  of  the  actual  non- 
linearity. This  gain  is  known  as  the  "describing  function".  For  a given 
probability  distribution  of  the  input,  (e.g.,  Gaussian,  as  shown  in  Figure  2), 
the  describing  function  is  determined  by  the  mean  and  standard  deviation  of 
the  input.  This  method  of  approximating  nonlinearities  is  known  as  quasi- 
linearization. 

To  apply  the  covariance  propagation  equations  1.4,  1.5  to  a nonlinear 
system  each  nonlinearity  is  replaced  by  its  describing  function,  which  is 
computed  for  the  mean  and  standard  deviation  of  the  input  at  the  beginning 
of  each  computing  interval.  By  updating  these  describing  functions  at  each 
computing  interval,  good  accuracy  can  be  obtained  even  for  complicated 
nonlinearities . 

We  shall  refer  to  this  technique  as  "nonlinear  covariance  propagation" 
(N.C.P.)  although  all  the  calculations  performed  within  any  computing 
interval  are  linear.  The  N.C.P.  technique  was  originally  developed 
for  analyses  of  missile  guidance  and  control  systems  (References  4, 

5 and  6).  The  computer  program  described  in  these  references  (the  CADET 
program)  was  not  used  as  the  basis  for  Vought's  launch  and  recovery  dynamics 


4.  Warren,  R.  S.  and  J.  Siegel,  SAM-D  Performance  Analyses  Usina  CADFT. 

Vols  I and  II.  TASC  TR  268-1,-2,  Contract  DAAH01-72-C-0826.'’  Aug.  73. 

5.  J.  H.  Taylor,  "Handbook  for  the  Direct  Statistical  Analysis  of  Missile 
Guidance  Systems  via  CADET,"  The  Analytic  Sciences  Corp.,  TR  385-2,  ONR 
Contract  No.  N00014-73-C-0213,  May  1975. 

6.  J.  H.  Taylor  and  C.  F.  Price,  "Direct  Statistical  Analysis  of  Missile 
Guidance  Systems  via  CADET,  Final  Report,"  ONR  Contract  No.  N00014-73-C-0213, 
March  1976. 
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A = SMALL  MEAN, 
LARGE  VARIANCE 

. B = LARGE  MEAN, 

SMALL  VARIANCE 


LOCAL 


IS  THE 


DESCRIBING  FUNCTION1 


FIGURE  2 ILLUSTRATION  OF  QUASI-LINEARIZATION 


program.  Instead  an  alternative  N.C.P.  program  'COVAR'  (Reference  7)  was 
selected  as  being  better  suited  for  our  technical  requirements.  The  N.C.P. 
program  presented  in  this  report  is  specially  written  for  modeling  the 
VSTOL  launch  and  recovery  system  and  incorporates  numerous  changes  from 
COVAR. 

These  changes  were  necessary  because,  compared  with  the  missile  guidance 
systems  studied  in  References  4,  5,  6,  and  7,  the  launch  and  recovery  system 
involves  more  complex  aerodynamic  nonlinearities,  a larger  number  of  state 
variables,  and  some  elements  which  require  special  modeling  techniques  (e.g., 
the  human  pilot).  Table  1 compares  the  CADET  program  versus  COVAR,  and 
VOLAR  (Vought  Launch  arid  Recovery  Dynamics  Program).  CADET,  COVAR,  and  VOLAR 
offer  alternative  approaches  for  solving  the  same  problem,  i.e.  evaluation 
of  a time  varying  function  P.  Main  advantages  of  VOLAR,  not  indicated  in 
Table  1,  are  convenient  input/output  features  and  a programming  philosophy 
that  facilitates  the  many  varied  models  required  for  launch  and  recovery 
analysis. 

1.5  Organization  of  this  Report 

This  report  is  arranged  so  that  each  section  leads  step-by-step  through 
two  example  calculations  of  VSTOL  launch  and  recovery  dynamics.  Both  examples 
relate  to  the  same  airplane  and  ship,  i.e.,  an  AV-8A  landing  on  an  LPD-class 
ship.  This  airplane-ship  combination  is  selected  because  flight  test  data 
are  available  to  check  the  predictions  of  the  VOLAR  Dropram. 

A ey  difference  between  the  two  examples  relates  to  the  human  pilot 
models  employed.  The  first  example  employs  a relatively  simple  pilot  model, 
the  parameters  of  which  are  determined  by  verbal  rules.  This  is  called  the 
"classical"  pilot  model.  In  the  second  example  the  pilot  model  parameters  are 
determined  automatically  as  the  result  of  a mathematical  optimization 
procedure.  This  model  is  known  as  the  "optimal  pilot"  model.  Both 
models  are  equivalent.  The  choice  of  model  is  a matter  of  user  preference  and 
convenience.  The  former  model  is  computationally  simple  but  demands  user  skill 
and  experience  to  select  valid  pilot  parameters.  The  optimal  controller  model 
avoids  this  restriction,  but  reouires  a considerable  increase  in  computer 
storage  and  run  time. 

% 

7.  Abzug,  M.  J.,  User's  Manual  for  FORTRAN  COVAR  Computer  Programs, 

Aeronautical  Consultant  Associates,  14951  Camarosa  Dr.,  Pacific  Palisades, 
Ca.  90272,  Rpt.  ACA  R-122.  Jan.  1976. 
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TABLE  1.  Comparison  of  Nonlinear  Covariance 
Propagation  Programs. 


CADET 


COVAR 


VOLAR 


ADVANTAGES/ 

DISADVANTAGES 


Treatment  of 
analytic 
nonlinearities, 
such  as  sin  e 


Treatment  of 
discontinuous 
nonlinearities , 
|such  as  limiters 
and  dead  zones 


Mean  solution 


Covariance 

Solution 


Gaussian 

describing 

functions 

(6DF) 


Linearization 

around 

reference 

trajectory 

found  using 

nonlinear 

equations 


Same  as 
COVAR 


COVAR  needs 
fewer  GDF's, 
reducing 
complexity 
of  program 


GDF 


Same 


Same 


Methods  are 
the  same 


Propagated 
linearly, 
using  GDF 
in 


Propagated 
nonli nearly, 
using 


Propagated 
nonli nearly 
using 


n 


Fm  + Gw  m = f (m)  + Gw  m = f(m)  + Gw 


t 


Coupled 


Input 


t 


COVAR  and 
VOLAR  provide 
more  accurate 
mean  solution; 
in  COVAR  mean 
solution  is 
offline  reduc 
ing  program 
size 


-Interacti  ve 


I 


I 


i 


Propagated  Propagated  Propagated 

linearly,  linearly,  using  linearly 

using  GDF  GDF  with  means  as  in  CADET 

in  found  in  non- 

T T linear  solution 
P=FP+PFi+GQGi  and 

P=FP+PFT+GQGT 


CADET  and 
VOLAR  have 
capability  to 
provide  more 
accurate 
covari ances 
for  large 
perturbations . 


NOTE:  m * mean  vector  P 

F = system  state  matrix  0 

f * system  nonlinear  state  equations 

G = noise  weighting  matrix  w 


covariance  matrix 
white  noise  power  spectral 
density  matrix 
white  noise  inputs 


♦VOLAR  may  be  'input'  or  'coupled'  to  suit  the  individual  requirement. 
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Section  2 of  this  report  presents  a aeneral  formulation  of  the  aircraft 
equations  of  motion  in  a form  compatible  with  Eauation  1.1,  known  as  "state- 
variable"  form.  For  some  launch  and  recovery  situations  the  lonoitudinal  and 
lateral  aircraft  degrees  of  freedom  may  be  effectively  coupled  (e.q.,  for 
aircraft  with  large  asymmetries  such  as  helicopters,  or  aircraft  with  high 
turn  rates  or  large  mean  sidesliD  angles).  Hence  Section  2 presents  the  fully 
coupled  algebraic  equations  of  motion  in  state  variable  form.  The  uncoupled 
equations  are  also  presented. 

Section  3 discusses  the  selection  of  AV-8A  data  and  presents  some  general 
information  on  modeling  stability  augmenter  systems,  usina  the  AV-8A  as  an 
example.  Section  3 concludes  with  a summary  of  the  AV-8A  numerical  data  used 
in  the  subsequent  example  calculations. 

Section  4 describes  a model  for  the  air  wake  of  a LPD-class  ship.  This 
model  is  based  upon  wind  tunnel  measurements  of  a FF1052  frigate  with  appropriate 
scaling  corrections  applied. 

Section  5 summarizes  the  classical  pilot  model.  Section  6 describes  how 
this  model,  plus  the  air  wake  and  AV-8A  models,  is  inputted  to  the  VOLAR 
program.  The  resulting  time  histories  of  means  and  variances  of  key  state 
variables  computed  by  VOLAR  are  presented  in  Section  6.  This  first  VOLAR 
example  applies  for  an  approach  performed  in  a low  sea  state,  for  which  ship 
motions  are  negligible. 

The  "optimal  controller"  pilot  model  is  described  in  detail  in  Section  7. 

In  Section  8 this  model  is  applied  to  compute  means  and  variances  of  the 
longitudinal  degrees  of  freedom  of  an  AV-8A  performing  an  approach  to  an  LPD- 
class  ship,  in  a low  sea  state  with  20  kts  wind  over  deck. 

Section  9 describes  the  ship  motion  model  that  Vought  has  developed  from 
the  Navy  RAOH  model.  In  Section  10,  the  example  of  Section  8 is  generalized 
to  include  ship  motion.  This  example  assumes  a moderate  sea  state  with  35  kts 
wind  over  deck. 

Section  11  summarizes  the  limited  data  currently  available  on  ground 
effects,  or  (more  properly)  'proximity1  effects.  Section  12  covers  some 
miscellaneous  topics  relating  to  covariance  propaaation  methods,  with  particular 
reference  to  interpretation  of  results  in  terms  of  probability  of  successful 
launch  and  recovery.  In  Section  13  the  results  of  the  example  calculations 
of  Sections  6,  8,  and  10  are  comDared  with  flight  test  data. 
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Conclusions  are  given  in  Section  14.  A User's  Manual,  program  listing, 
and  some  specialized  mathematical  topics  are  presented  in  appendices. 


2.0  AIRCRAFT  DYNAMIC  MODEL 
AIRCRAFT  EQUATIONS  OF  MOTION  IN  STATE-VARIABLE  FORM 

2.1  Introduction 

The  launch  and  recovery  dynamics  program,  like  other  computational  methods 
based  on  modern  control  theory,  requires  the  rystem  dynamic  equations  to  be 
written  in  the  form 

lx  = Fx  + Gw  (2.1) 


where 

x = the  system  output  state  vector 
F = the  system  matrix 
w = the  input  state  vector 
G = the  input  matrix 
1 = the  identity  matrix 

The  above  form  will  be  referred  to  here  as  the  "state-variable 
form".  It  differs  from  the  form  in  which  airplane  equations  of  motion  are 
conventionally  written.  The  most  common  form  for  airplane  equations  involves 
second  derivatives  of  an  output  state  vector  x',  thus 


d2x* 

dt2 


dx' 

dt 


Cx' 


Dw' 


(2.2) 


where 

w'  = control  and  gust  input  vector 

By  introducing  additional  state  variables  augmenting  x'  to  x"  this  can  be 
reduced  to  a first  order  equation  of  the  form 


M 


dx" 

dt 


+ 


Nx" 


Ew' 


(2.3) 


In  theory,  reduction  to  the  state-variable  form  is  possible  by 
substituting  numerical  values  in  equation  2.3  and  subsequently  premultiplying 
by  M“l,  however,  it  is  generally  more  convenient  and  more  accurate  to 
re-formulate  the  original  literal  equations  in  the  algebraic  form  of 
equation  2.1  before  substituting  numerical  values.  This  section  describes 
the  appropriate  literal  forms  for  the  longitudinal  and  lateral 
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equations  of  motion,  and  discusses  the  reduction  of  cross-coupled  lateral- 
plus-longitudinal  equations  of  motion  to  state-variable  forms. 

2.2  Transformation  of  Lateral  Equations  to  State-Variable  Form 


From  Reference  8 the  linearized  equations  describing  small 
lateral  perturbations  from  straight-line  flight  are: 


mY,IxL,IyN,  are  the  applied  aerodynamic  forces  and  moments  due  to  motion 
perturbations  v,  p,  r and  also  due  to  control  deflections.  It  is  generally 
assumed  that  the  aerodynamic  forces  and  moments  due  to  v,  p,  r,  and  higher 
order  derivatives  are  negligible.  With  this  assumption,  equation  2.4  can 
be  reduced  to  the  desired  form  through  the  following  manipulations:  multiply 
the  third  row  by  Ixz/I  and  add  the  result  to  the  second  row;  multiply  the 
second  row  by  Ixz/Iz  and  add  the  result  to  the  third  row;  multiply  the  fifth 
row  by  sin  0 and  add  to  the  fourth  row;  finally,  divide  the  second  and 

0 I 2 

third  rows  by  (1  - xz  ).  These  manipulations  yield: 

yz 

8.  McRuer , P.  T.,  T.  |_.  Ashkenas,  and  F.  D.  Graham,  Aircraft  Dynamics  ard 
Automatic  Control,  Princeton  University  Press,  Princeton,  New  Jersey, 
1973. 
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where  L ' , N'  are  prime  derivatives  as  defined  in  Reference  8: 
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Li  + Ni  T 


xz 
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1 - 


■y r 


v - 
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Vz 


(2.6) 


With  the  addition  of  the  kinematic  relation  for  lateral  displacement  from  the 
desired  flight  path,  Ay,  equation  (2.5)  can  be  written  in  the  desired  form 


[F] 


[I]  (x) 

[F]  {x> 

+ [G]  (w) 

(2.7) 

(X  }T  = 

CL 

* 

> 

r,  <f> , ip.  Ay} 

= 

\ 

W+Y 

0 p 

-U  +Y  g 

or  s 

COS 

0o 

0 

0 

L ' 

L ' 

L ' 

0 

0 

0 

V 

P 

r 

N ' 

N ' 

0 

0 

0 

(2.8) 

V 

P 

r 

0 

1 

tan  0Q 

0 

0 

0 

0 

0 

1 

0 

0 

0 

COS  6 „ 
n 

Un  cos 

0O 

_ 1 

0 

0 

0 

+w0  sin 

0 

0o 

0 _ 

14 


[G] 


{w}1 


Y6A 

Y6R 

-Yv 

0 

0 

0 

L6a' 

L6r‘ 

-v 

0 

0 

0 

V 

V 

-v 

0 

0 

0 

0 

0 

0 

0 

o' 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

V 

•* 

o 

CD 

> 

0} 

(2.9) 


(2.10) 


This  formulation  excludes  the  aerodynamic  forces  due  to  lateral  gust  gradient 
terms,  9v  , 9u  ; these  are  expected  to  be  negligible  for  VSTOL  aircraft  with 

■ J — — i M 

9 9 

x y 

closely  grouped  engines  at  low  speeds. 

2.3  Longitudinal  Equations  in  State-Variable  Form' 

From  reference  8 the  linearized  equations  describing  small 
perturbations  from  straight-line  flight  are: 
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(2.11) 


where  mX,  mZ,  IyM  are  the  applied  aerodynamic  forces  and  moments  due  to  the 
controls  and  gust  inputs,  and  to  the  motion  perturbations.  It  is  generally 
assumed  that  X,  Z,  and  M are  .linearly  dependent  on  u,  w,  q,  and  also  on  the 
time  derivative  of  w,  which  is  assumed  to  affect  the  aircraft  response  via 
derivatives  X^,  Z^,  M^.  Usually  only  is  retained,  yielding 


(2.12) 
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where  M = M - M^.  w = Mqq  + l^w  + M^u 

To  reduce  equation  (2.13)  to  the  desired  form  multiply  the 
second  row  by  M^,  and  subtract  the  result  from  the  fourth  row.  To  this  is  adjoined 
kinematic  equations  for  the  horizontal  and  vertical  position  increments  Ax,  Az. 

This  yields  an  equation  of  the  form  of  equation  (2.1). 

Ix  = Fx  + Gw  (2.14) 

where  {x}  ^ = {u,  v,  e,  q.  Ax,  Az) 
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6^ , 62,  63  represent  the  deflections  of  the  appropriate  longitudinal 
controls,  e.g.  stick,  throttle  vectoring  nozzles,  etc. 
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2.4  Transformation  of  Cross-Coupled  Equations  of  Motion  to 

State-Variable  Form 

The  longitudinal  and  lateral  equations  of  motion  become  coupled  for 
aircraft  performing  small  perturbations  from  curved  flight  paths  and  also 
for  straight-line  flight  of  aircraft  with  asymmetric  aerodynamic  and/or 
inertial  characteristics,  e.g.  helicopters.  The  transformation  to  state- 
variable  form  of  the  standard  coupled  equations  of  motion  (e.g.,  as  given  in 
Reference  8)  is  more  complicated  than  for  the  uncoupled  cases  considered 
previously.  The  procedure  described  below  requires  that  the  aircraft  x and  y 
body  axes  shall  be  defined  such  that  I = I = C,  but  I is  in  general  non- 

J\j  j L XZ 

zero.  The  assumption  is  made  that  all  acceleration  derivatives  such  as 
can  be  neglected. 

For  compactness  the  eauations  will  be  presented  here  in  terms  of  the 
uncoupled  longitudinal  and  lateral  equations  already  derived,  with  additional 
terms  added  as  necessary. 

The  12-component  state  vector  {x}  is  defined  as 


{x}  = (u,  w,  0,  q,  Ax,  Az;  v,  p,  r,  $,  ip.  Ay} 

The  (n  + 3)  component  input  vector  is  defined  as 


(2.17) 


*ug’  Wg’  6l  ^2  • " • • 6p;  vg»  6p+i  •••  6n} 


(2.18) 


where  u , v , and  w are  gust  inertial  velocities,  6j,  62  ...  6 represent 
y y y 

longitudinal  control  deflections,  and  6p+-j , 6p+2 » • • • 6n  represent  lateral 
control  deflections. 

For  these  state  and  control  vectors,  the  equations  of  motion  in  the 
desired  form  are 

lx  = Fx  + Gw  (2.19) 

where  F and  G are  defined  as  described  below. 

DEFINITION  OF  THE  F-MATRIX 

The  F-matrix  is  a 12  x 12  matrix  with  rows  compatible  with  the  x vector 
defined  previously. 
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The  order  selected  for  the  elements  listed  in  Equation  2.17  facilitates 
decomposition  of  the  F-matrix  into  the  sub-matrices  Fj^,  F^.  f^l*  ^22  as 
shown  below,  which  for  P0  = RQ  = Qo  and  zero  values  of  cross-coupling 
derivatives  such  as  Lq,  Mr,  etc.  yields  the  uncoupled  lateral  and  longitudinal 
F-matrices  previously  discussed  as  Fjj,  F£2,  with  F^2  = F21  = 


LONG  j 

f 

f Fn  ! 

F12 

jxL0NG< 

! 

i + 

Gn  ! 

! 

LAT  ‘ 

! 

F21  j 

F22  _ 

LXLAT  - 

r T 

i 

_ G21  J 

The  size  of  the  F-matrix  precludes  single-page  presentation.  Therefore, 
to  facilitate  the  visualization  of  physical  effects  the  F-matrix  is  given  in 
transposed  form  in  Table  2,  so  that  all  the  contributions  to  any  particular 
given  force  or  moment  appear  on  a single  page.  For  compactness  the  abbrevia- 
tions s = sin,  c = cos,  t = tan  are  employed,  and  roll  and  yaw  accelerations 
are  denoted  by  L ' , L",  N',  N",  as  explained  telow.  The  G-matrix  is  given  in 
Table  2. 

DEFINITIONS  OF  L ' , N',  L",  N"  TERMS 


I. 


L.+fi 


xz 


' ' >* 


Vi 


XZ 


1 - 


XZ  ' 


Ni' 


<[*z> 


where  i = u,w,v,p,q,  r 

I-  I.-I,  . r »„)2, 

% ♦ <Np ♦ -V*  V -fj  / D - T-ff~] 


L " = Tl 

P L c 


xz 


x z 
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TRANSPOSED  F-MATRIX  (CONTINUED) 


3.0  AIRCRAFT  DATA  AND  STABILITY  AUGMENTER  SYSTEM  MODEL 


3 .1  Introduction 

It  is  desirable  to  demonstrate  the  validity  of  the  nonlinear  covariance 
propagation  technique  by  comparing  the  results  that  it  yields  versus  data 
from  actual  flight  operations.  This  requires  mathematical  models  of  the  ship 
motion,  the  airwake  and  the  atmospheric  conditions  prevailing  during  the 
flight  tests,  plus  a model  of  the  airframe  and  control  system  dynamics.  The 
data  currently  available  do  not  permit  one-to-one  matching  of  all  these 
factors.  For  example,  published  quantitative  flight  data  on  the  AV-8  air- 
craft are  extremely  limited,  and  not  compatible  with  the  available  ship 
motion  and  airwake  models.  Thus  dynamic  interface  charts  are  available  for 
the  LPD  (Reference  9)  and  LPH  class  ships  (Reference  10),  but  satisfactory 
airwake  models  based  on  wind  tunnel  data  are  only  available  for  the  FF-1052 
and  DD-963  ships  (References  11  and  12  respectively).  It  is  therefore  not 
advisable  to  focus  narrowly  on  a super-accurate  model  of  the  airframe  while 
tolerating  gross  uncertainties  in  other  system  elements.  This  is  particularly 
true  for  simulations  which  replace  the  human  pilot  by  a mathematical  model. 

The  human  pilot  may  require  modeling  of  secondary  dynamic  effects  (such  as 
stall  and  buffet  frequency  characteristics)  in  order  to  reassure  him  of  the 
overall  fidelity  of  the  simulation,  but  such  detail  is  unnecessary  if  a pilot 
model  is  employed.  These  considerations  influence  the  selection  of  aircraft 
data  from  the  sources  listed  below. 
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AV-8  aircraft  models  of  sufficient  detail  for  correlation  have  been 
published  in  References  13,  14,  and  15.  All  these  models  are  based  on  McAir 
data.  The  model  of  Reference  13  was  developed  for  the  AV-8A  and  includes 
numerous  nonlinearities  which  are  represented  by  tables  of  the  appropriate 
dependent  variable  for  given  values  of  one  or  at  most  two  independent 
variables.  Graphs  of  these  tables  are  included  in  Reference  13,  the  tabulated 
values  of  the  independent  variables  being  joined  by  straight  lines.  Reference 
14  presents  the  basic  McAir  aerodynamic  data  on  the  AV-8A  and  AV-8B  and 
Reference  15  presents  AV-8B  models  used  for  the  X-22A  in-flight  simulator. 

Data  on  the  AV-8A  stability  augmenter  system  (S.A.S.)  is  given  in  Reference 
14.  S.A.S.  was  employed  during  the  flight  tests  described  in  References  9 
and  10.  After  comparing  the  available  data  sources,  Reference  15  was 
selected  as  the  most  suitable  for  our  use,  based  on  its  demonstrated  accuracy, 
and  its  simplicity  compared  to  the  alternative  models.  The  model  represents 
the  AV-8B  rather  than  the  AV-8A,  but  for  the  flight  conditions  of  interest 
here,  the  differences  are  minor,  as  noted  in  Reference  14,  and  can  be 
corrected  as  indicated  in  the  following  subsection. 

3.2  AV-8A  and  AV-8B  Models 

Reference  15  lists  aerodynamic  data  for  two  AV-8B  approach  profiles,  one 
for  a 30  degree  approach  from  105  knots,  the  other  for  a 5 degree  approach  from 
65  knots  and  also  gives  "single-model"  data  which  approximate  the  aircraft 
characteristics  for  both  profiles.  Table  4 lists  these  "single-model  data", 
which  are  used  in  the  present  simulation.  Although  some  nonlinear  effects 
have  been  neglected,  these  derivatives  are  shown  in  Reference  15  to  give 
good  agreement  with  the  results  of  the  McAir  fully  nonlinear  simulation. 


13.  Nave,  R.  L.,  Progress  Toward  a Computerized  VSTOL/Small  Platform  Landing 

Dynamics  Investigation  Model,  NADC  Report  77024-30,  Naval  Air  Dev.  Ctr., 
Warminster,  Pa.,  1977. 

14.  Anon.,  AV-8B  Simulator  Evaluation,  VST0L  Performance,  McDonnell  Douglas 

Corp.  Rept.  MDC-A3922,  Contract  N00019-75-C-0487,  March  1976. 

15.  Lebacqz,  J.  V.,  Summary  Documentation  of  AV-8A  Model  Development  and 

X-22A  Simulation  of  AV-8B,  Calspan  Corp.,  Buffalo,  N.  Y. 

X-22A  TM  No.  98,  W/A  P63-054,  July  1977. 
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TABLE  4(a).  "SINGLE-MODEL"  AV-8B  LONGITUDINAL  MODEL 


VKT> 

0 

30 

50 

65 

80 

105 

Xu 

-.044 

-.044 

-.044 

-.044 

-.044 

-.044 

Zu 

0.0 

-.023 

-.054 

-.092 

-.101 

-.105 

Mu 

0.0 

-.0009 

-.0020 

-.0026 

-.0008 

0.0 

Xw 

0.0 

0.0 

+.0035 

+ .010 

+.0175 

+ .035 

Zw 

-.018 

-.125 

-.195 

-.240 

-.300 

o 

a> 

CO 

i 

Mw 

+.0042 

+ .0047 

+.0040 

+.0021 

-.0015 

-.0070 

Xq 

0 

+ .03 

+ .05 

+ .052 

+ .085 

.11 

Zq 

-.05 

-.29 

-.40 

-.47 

-.61 

-.79 

"o 

-.056 

-.118 

-.160 

-.191 

-.224 

-.280 

X ^ 
A6ES, 

-.161 

-.151 

-.145 

-.139 

-.134 

-.126 

Z5ES 

-.350 

-.340 

-.335 

-.330 

-.295 

-.205 

M6ES+ 

.230 

.235 

.237 

.240 

.241 

.243 

X6T 

.34 

.30 

.262 

.225 

.180 

.140 

V 

-2.55 

-2.46 

-2.37 

-2.26 

-2.10 

-1.75 

M6T* 

-0.036’ 

-0.036 

-0.036 

-0.036 

-0.036 

-0.036 

w *+ 

A6T 

.66 

.60 

Z6T** 

- 

- 

- 

-2.11 

- 

-1.41 

X * 

V 

-.555 

-.546 

-.530 

-.510 

-.476 

-.400 

z . 

V 

-.066 

-.061 

-.050 

-.028 

-.023 

o 

o 

1 

V 

0 

0 

0 

0 

0 

0 

★★ 

X3j 

- 

- 

- 

-.456 

- 

-.380 

y ★★ 

ej 

- 

- 

-.066 

-.055 

*9j  « 81  degrees  **For  trimmed  condition 

t Units  of  6ES  are  not  defined  in  Reference  15 


n 
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TABLE  4(b): 

"SINGLE- 

■MODEL"  AV-8B 

LATFRAL- 

-DIRECTIONAL 

MODEL 

IJ  I 

D 

• a 1 

Vo 

0 

30 

50 

65 

80 

105 

li  I 

Y 

-.063 

-.088 

-.104 

-.120 

-.138 

n 1 

V 

•»  1 

V 

Bum|I 

-.0144 

-.0197 

-.0204 

-.0184 

-.0104 

N ' 

-.0021 

-.0010 

0.0 

+.0014 

.0047 

n 1 

V 

11 

Yp 

0 

0 

0 

0 

0 

0 

n \ 

V 

-.13 

-.42 

-.62 

-.79 

-1.0 

-■  .38 

V 

-.005 

-.032 

-.053 

-.072 

-.091 

-.126 

Yr 

-.225 

- .24 

-.24 

-.235 

-.235 

-.215 

Lr' 

.015 

.15 

.24 

.31 

.375 

.49 

V 

-.042 

-.088 

-.12 

-.142 

-.164 

-.203 

Y6as  + 

0 

-.015 

-.026 

-.03? 

-.065 

Li  1] 

L6as't 

0.5 

0.5 

0.5 

0.5 

0.5 

n 

N5as'+ 

.030 

.030 

.031 

.033 

.035 

.040 

»*  1 

Y6RP  + 

-.68 

-.67 

-.70 

-.75 

-.80 

-.90 

• • ■ 

L6RP'+ 

-.065 

-.075 

-.080 

-.095 

-.13 

n I 

N6RP'+ 

.225 

.235 

.243 

.248 

.255 

.265 

11  I 

_ II 
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In  the  present  simulation,  nonlinearities  are  introduced  to  model  the 
limits  on  maximum  roll  and  yaw  acceleration  capabilities.  Reference  15 
notes  the  Nv'  derivative  as  being  very  nonlinear  at  certain  speeds.  The 
N'  versus  v variation  can  be  modeled  as  an  odd  function  of  v comprised  of 
linear  segments  with  break  points  at  -12,  -8,  8,  and  12  degrees  (see 
Figure  3).  This  was  not  required  because  Ny'  was  linear  over  the  speed 
range  simulated  here. 

AV-8A  maximum  control  power  capabilities  are  obtained  from  Reference 
14  as: 

lMAX  = rad/sec2,  MMAX  = 0.25  rad/sec2,  N^x  = 0.125  rad/sec2 


u 

o 


u 

n 


Due  to  the  bleed  system  the  maximum  T/W  ratio  depends  on  the  utilization  of 
the  above  control  moments.  The  appropriate  T/W  describing  function  could  be 
programmed  to  reflect  this  dependence,  but  this  does  not  appear  to  be 
necessary  for  purposes  of  demonstrating  correlation  with  the  particular  flight 
data  of  interest  here.  A mean  value  of  T/W  = 1.3  was  used,  appropriate  to  an 
installed  static  thrust  of  21,500  pounds  out  of  ground  effect. 

Reference  14  discusses  the  differences  in  stability  and  control  character- 
istics of  the  AV-8A  and  AV-8B.  For  the  low  speeds  of  interest  here  (<65  kt) 
the  effects  of  increased  span  and  area  of  the  AV-8B  are  apparently  secondary, 
compared  to  the  improvement  in  roll  control  of  the  AV-8B  due  to  realignment 
of  the  roll  puffers  to  reduce  adverse  yaw.  This  could  be  simulated  by 
adjusting  the  ratio  of  N6AS‘  to  if  desired.  A second  difference  involves 

ground  effects,  which  Reference  14  notes  are  improved  on  the  AV-8C5  due  to 
the  fuselage  strakes.  The  numerical  simulation  examples  presented  in  this 
report  are  mainly  out  of  ground  effect,  so  the  AV-8B  representation  is 
adequate.  Finally,  the  AV-8B  is  claimed  to  have  improved  inlet  stall 
characteristics  due  to  the  use  of  two  rows  of  blow-in  doors  near  the  inlet 
lips.  From  References  14,  15,  and  16  it  appears  that  this  change  does  not 
affect  L v ' except  at  large  sideslip  angles,  since  the  Lv'  data  supplied  in 
these  references  are  essentially  linear  for  moderate  sideslip  angles. 


16.  Lacy,  T.  R.,  MIL-F-83300;  View  from  an  Aircraft  Designer,  Proc.  Navy/NASA 
VSTOL  Flying  Qualities  Workshop,  April  1977,  Naval  Postgraduate 
School,  Monterey,  California,  Published  August  1977. 
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3.3  Representation  of  Stability  Auamenter  Systems 

The  representation  of  stability  augmenter  systems  for  the  Launch  and 
Recovery  Dynamics  Program  and  for  covariance  propagation  methods  in  general 
is  complicated  by  the  need  to  also  represent  the  human  pilot.  The  output  of 
the  human  pilot  model  must  be  summed  with  the  output  of  the  S.A.S.  in  order 
to  compute  the  total  forces  and  moments  applied  to  the  aircraft  by  the  controls. 
Complications  arise  because  the  authority  limits  on  the  S.A.S.  are  generally 
set  such  that  the  forces  and  moments  produced  as  a result  of  S.A.S.  commands 
are  considerably  less  than  the  maximum  values  of  the  forces  and  moments  that 
can  be  produced  in  response  to  pilot  commands. 

Figure  4 (page  32)  illustrates  these  considerations.  A fly-by-wire  system  is 
selected  for  simplicity  because  only  one  set  of  actuators  need  be  considered, 

(We  assume  that  pilot  commands  and  S.A.S.  commands  generate  electrical  signals 
which  are  summed  and  input  to  these  actuators).  The  principles  discussed 
below  apply  to  other  types  of  flight  control  systems,  but  their  implementation 
changes  for  systems  where  some  (but  not  all)  of  the  commanded  control  deflec- 
tions are  fed  through  some  or  all  of  the  actuators.  The  describing  functions 
shown  on  Figure  4 are  computed  as  follows. 

The  maximum  values  of  the  forces  and  moments  that  can  be  generated  by 
the  flight  control  system  determine  the  lower  describing  function  block  shown 
on  Figure  4.  For  VSTOL  aircraft  this  block  typically  represents  four 
'saturation'  describing  functions,  one  each  for  yaw,  pitch,  and  roll  and 
maximum  thrust.  The  upper  describing  function  block  represents  the  saturation 
describing  functions  that  define  the  limits  on  the  maximum  values  of  the  S.A.S. 
commands.  For  covariance  propagation  it  is  necessary  to  sum  the  output  of  the 
pilot  model  (upi-lot)  with  the  linearized  output  of  the  S.A.S.  describing 
function  (u$).  The  mean,  Uj,  and  variance,  °uj,  of  the  sunrned  signal, 
uj,  determine  the  effective  gain  of  the  lower  describing  function  in 
Figure  4. 

3.4  AV-8A  Autopilot 

Later  in  this  report  the  Launch  and  Recovery  Dynamics  Program  is  employed 
to  compute  the  behavior  of  the  AV-8A  approaching  and  landing  on  a small 
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ship  of  the  LPD  class.  The  program's  predictions  are  correlated  with  flight 
test  results  for  the  AV-8A  landing  on  the  U.S.S.  'Raleigh'.  These  flight  test 
data  were  obtained  on  early  AV-8A's  which  were  fitted  with  a roll  S.A.S.  of 
very  limited  authority  and  had  no  pitch  or  yaw  stability  augmentation.  On 
subsequent  MSA's  the  S.A.S.  was  upgraded  to  a level  similar  to  that  of  the 
AV-8B  S.A.S.  described  in  Reference  14  except  that  no  roll-yaw  interconnect 
was  fitted.  A brief  summary  of  this  S.A.S.  is  given  below.  This  summary 
facilitates  the  extrapolation  of  the  results  given  in  this  report  to  later 
versions  of  the  AV-8A  and  also  to  the  AV-8B. 

Block  diagrams  of  the  AV-8A  S.A.S:.  are  given  in  Reference  14.  These  may 
be  expressed  in  the  form  of  Figure  4 with  the  following  substitutions. 


TABLE  5.  AV-8A&B  S.A.S.  DYNAMICS 


STATE 

VARIABLE 


ACTUATOR 

TIME  CONSTANT  (SEC 


S.A.S.  COMPENSATION 


deg  6 


K = 0.669 


0.0175 


deg/sec 


deg  6 


K = 0.83 


0.0175 


ft/sec 
deg  6a 


K = 0.346 


deg  6 


0.01755 

(front  reaction 
control  system) 
0.0833 

(aft  reaction 
control  system 
plus  stabilizer) 


K = 0.495 


TABLE  6.  AV-8A&B  S.A.S.  CONTROL  AUTHORITY 


CONTROL 


UT 
LIMITS 


±1.5  Deg.  Stabiliator 
±2.0  Deg.  Ailerons 
±5.0  Deg.  Rudder 


0.266 


0.123 

0.50 


For  details  of  the  roll -yaw  interconnect  system  which  is  fitted  to  the 
AV-8B  only  see  Ref.  14.  The  very  low  authority  of  the  roll  S.A.S.  should 
be  noted. 


V °us 
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Figure  4.  Simplified  Representation  of  S.A.S.  for  Fly-By-Wire 
System 


4.0  ATMOSPHERIC  MODEL 


4.1  Introduction 

The  Vought  launch  and  recovery  dynamics  program  is  compatible  with 
many  different  models  of  the  atmospheric  environment.  The  simplest  type 
of  model  applies  when  the  airplane  is  completely  outside  the  ship  airwake. 

For  this  case  the  atmospheric  model  comprises  steady  wind  (including  wind 
shear),  and  turbulence.  For  critical  phases  of  VSTOL  operations  the 
airplane  is  either  within  or  near  the  physical  airwake;  this  necessitates 
a more  complicated  model,  such  as  that  presented  below.  This  model  is, 
however,  not  restricted  to  the  space  within  the  physical  airwake,  since  it 
blends  into  a general  low  altitude  atmospheric  model  as  the  airwake 
boundaries  are  crossed. 

The  ship  airwake  model  is  derived  from  wind  tunnel  measurements  on  a 
model  of  an  FF  1052  frigate.  It  was  previously  presented  in  Reference  11, 
so  only  a brief  description  is  given  here.  The  model  gives  airwake 
properties  as  functions  of  the  airplane  position  with  respect  to  the 
ship,  for  six  wind  over  deck  conditions.  These  six  conditions  are  presented 
in  Table  7.  VOLAR  interpolates  between  these  six  conditions  to  provide 
ship  airwake  information,  at  any  point  in  space,  for  a wind  over  deck 
orientation  between  30  degrees  and  50  degrees  to  port,  with  a magnitude 
between  20  knots  and  50  knots. 

The  major  departure  from  Reference  11  occurs  in  the  bookkeeping  of 
the  airwake  and  ambient  winds.  VOLAR  allows  the  user  to  select  no  airwake 
at  all  or  to  model  an  airplane  flying  in  and/or  out  of  the  airwake. 

Hence,  VOLAR  bookkeeping  is  different  from  that  of  Reference  11.  In  VOLAR, 
subroutine  AIRJAKE  returns  tne  incremental  disturbance  velocities  due  to 
the  ship's  presence,  not  the  total  airwake  velocities  (which  include  the 
ambient  wind)  given  by  Reference  11.  In  other  words,  the  ambient  wind  is 
subtracted  from  the  airwake  and  accounted  for  separately. 

VOLAR  can  also  represent  Dryden  type  turbulence  alone  or  in  addition 
to  the  ship  airwake  turbulence.*  In  the  work  reported  here,  Dryden  turbulence 
was  used  in  conjunction  with  the  ship  airwake  to  model  random  atmospheric 
disturbances  to  the  airframe. 

*V0LAR  is  not  restricted  to  Dryden  type  turbulence.  Any  turbulence  that 
can  be  represented  by  white  noise  fed  through  a physically  realizable  filter 
with  a finite  number  of  poles  and  zeros  can  be  modeled  by  VOLAR. 
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Table  7.  Wind-Over-Deck  Conditions  for 
Airwake  Model 

n 

ii 


VW0D  * KTS 

*s  -v  DEGREES 

11 

20 

-30 

20 

-50 

11 

35 

-30 

35 

-50 

ii 

45 

-30 

45 

-50 

l i 

is  the  wind  over 

deck  direction  with 

respect  to  the  ship  heading,  positive 

if  wind  over  deck 

is  coming  into  the 

starboard  side. 

4.2  Summation  of  Airwake  Components 

r B 

This  section  outlines  tiie  development  of  the  aircraft  airspeed,  V^, 
from  the  aircraft  body  axis  inertial  velocities  and  the  various  atmospheric 
wind  and/or  turbulence  models. 

4.2.1  Ship  Airwake  Model 

A schematic  diagram  depicting  the  ship  airwake  model  is  presented  in 
Figure  5.  This  flow  chart  shows  the  inputs,  summarizes  the  manipulations, 
and  shows  the  outputs  of  the  ship  airwake  calculations  in  VOLAR. 

n 
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Figure  5.  Airwake  Model  Flow  Chart 


4.2.2  Dryden  Turbulence  Model 

As  shown  in  Figure  6 » VOLAR  models  the  Dryden  turbulence  components  as 
outputs  of  first  order  filters  driven  by  white  noise.  One  filter  is  required 
for  each  axis,  and  the  noise  sources  are  uncorrelated. 

4.2.3  Wind  Velocities  at  Aircraft 

The  build-up  of  the  total  wind  velocity  is  depicted  in  Figure  7. 

From  the  aircraft  body  axis  components  of  the  wind  (Uwind,  Vwinc(,  Wwind)  the 
airspeed  components  are  readily  determined,  as  follows: 


Uas  “ U Uwind 
Vas  = V • Vwind 
Was  - W = “wind 


and  V.  = U2  + V2  + W2 
A as  as  as 


Independent 
White  Noise 
Sources 


Airplane  Body  Axis 
Random  Disturbances 


FIGURE  6.  TURBULENCE  SIMULATION 
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Figure  7.  Summation  of  Airwake  Components 
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5.0  CLASSICAL  PILOT  MODEL 
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5.1  Types  of  Pilot  Model 

The  human  pilot  is  an  essential  element  of  the  VSTOL  launch  and  recovery 
system  and  his  dynamics  should  be  modeled  as  accurately  as  anv  other  major 
element  of  the  system.  Two  aoproaches  to  modeling  the  human  Dilot  are 
currently  in  widesDread  use.  These  will  be  desionated  the  "classical"  model 
and  the  "optimal  pilot"  model  respectively.  The  dynamic  models  that 
result  from  these  two  modeling  methods  are  eoui valent  (if  both  models  are 
accurate) . 

The  rationale  for  choosino  one  pilot  model  over  the  other  will  become 
clear  after  reading  this  section  and  section  7.  These  sections 
describe  both  models. 

5.2  Example  of  Classical  Pilot  Model 

Since  the  classical  pilot  model  for  closed-loop  control  has  been 
documented  extensively  in  Reference  17  and  elsewhere  onlv  a highly  condensed 
description  will  be  given  here.  The  model  essentially  consists  of  blocks  with 
describing  functions  of  the  form 

Output  _ -ts  . <TLS  * ^ 

Input  " p (Tns  + 1) 

where 

Kp  = pilot  gain,  t = time  delav,  T^  = lead  time  constant, 

Tj^j  = neuromuscular  time  constant 

Additional  leads  and  lags  mav  be  introduced  as  necessarv.  These  leads  and 
lags  and  the  pilot  gain  are  adjusted  to  meet  certain  "verbal  adjustment  rules" 
given  in  Reference  17.  These  rules  are  unfortunately  not  stated  in  forms 
suitable  for  computer  mechanization.  Therefore  use  of  the  classical  pilot 
model  reauires  judgement  and  experience,  i.e.,  it  is  something  of  an  "art". 
This  comment  also  applies  to  the  choice  of  Iood  structure  into  which  the 


0 

n 


17.  McRuer,  D.  T.,  and  H.  R.  Jex,  A Review  of  fluasi-Linear  Pilot  Models, 
Trans.  l.E.E.E,  Vol . HFE-8,  No.  3,  SeDtember  1967 . 


38 


pilot  describino  function  blocks  are  inserted  for  multilooD  control  situations. 
Further  adjustments  are  reauired  to  correct  the  Dilot  model  for  variations  in 
system  dynamics  caused  bv  changes  in  f 1 i ah t condition.  For  lateral  control 
of  the  AV-8A  at  low  SDeeds  and  hover  the  variations  are  sufficiently  small 
that  a constant  loop  structure  and  fixed  oains  and  time  constants  may  be  used, 
as  shown  in  Figure  8.  The  numerical  values  shown  in  Figure  8 are  selected 
to  match  the  classical  oilot  model  verbal  adjustment  rules. 

Details  of  classical  pilot  models  for  other  VSTOL  control  situations 
(e.g.,  longitudinal  control)  are  oiven  in  References  18,  19  and  20. 


18.  Craig,  S.  J.,  and  A.  Campbell,  "Analysis  of  VTOL  Handling  Dualities 

Reauirements"  Part  I:  "Lonaitudinal  Hover  and  Transition,"  AFFDL- 
TR-67-179,  1968. 

19.  Craig,  S.  J.,  A.  Campbell,  and  R.  H.  Klein,  "Analysis  of  VTOL  Handling 

Dualities  Reauirements,  Part  II:  Lateral -Directional  Hover  and 
Transition,"  AFFDL-TR-67-179,  Part  II,  1970. 

20.  Craig,  S.  J.,  I.  L.  Ashkenas,  and  R.  K.  Heffley,  "Pilot  Background  and 

Vehicle  Parameters  Povernino  Control  Techniaues  in  STOL  Approach 
Situations,"  FAA  RD  72-69,  June  1972. 
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6.0  EXAMPLE  APPROACH  FOR  AV-8A  WITH  CLASSICAL  PILOT  MODEL 


6.1  Introduction 

To  familiarize  the  reader  with  the  VOLAR  program  as  soon  as  possible,  we 
now  present  a simple  example  of  its  use.  This  example  models  the  complete 
launch  and  recovery  system  by  combining  the  analytic  models  for  the  AV-8A, 
airwake,  and  human  pilot  described  in  Sections  3,  4,  and  5.  To  illustrate 
the  effects  of  the  limited  maximum  available  control  moments  the  example  is 
run  twice,  once  with  and  once  without  these  limits.  The  importance  of  these 
nonlinear  elements  can  be  assessed  by  comparing  the  two  sets  of  results. 

The  approach  flight  condition  selected  here  is  a low  speed  approach  at 
constant  50  feet  altitude  directly  into  a 35  kt  wind-over-deck,  as  shown  in 
Figure  9.  The  aircraft  initial  speed  relative  to  earth  xg  (see  Figure  9) 
is  66.516  fps  = 39.4  kt  = 20.28  m/s.  During  the  first  16  seconds  of  the 
simulated  approach  the  aircraft  decelerates  at  a constant  rate  until  it  is 
stationkeeping  with  the  ship.  The  VOLAR  simulation  covers  20  seconds  of 
flight,  the  last  four  seconds  representing  stationkeeping.  The  altitude 
(50  ft)  is  constant  throughout;  the  final  descent  stage  of  the  approach  is 
not  simulated  in  this  example,  although  it  is  simulated  in  examples  presented 
later  in  this  report.  A number  of  assumptions  are  made  for  this  particular 
example  which  simplify  the  calculations,  at  the  cost  of  restricting  the  range 
of  applicability  of  the  results.  These  assumptions  include: 

(1)  Ship  motion  effects  are  neglected. 

(2)  The  mean  sideslip  velocity  of  the  airplane  is  small. 

This  is  in  accord  with  AV-8A  operating  procedure.  It 
also  permits  the  aircraft  equations  of  motion  to  be 
separated  into  longitudinal  and  lateral  sets  of 
equations  which  can  be  solved  independently.  For  the 
present  example  only  the  lateral  set  is  considered,  the 
longitudinal  equations  are  solved  in  examples 
presented  in  Sections  8 and  10. 

The  aircraft  inertial  flight  path  is  curved.  This  follows  from 
Figure  9,  because  the  aircraft's  component  of  inertial  velocity  parallel 
to  the  xr-axis  decreases  with  time,  but  the  ship  inertial  lateral 
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FIG.  9 COMMANDED  LATERAL  APPROACH  PROFILE 


velocity  is  constant  = 20  cos  30°  = 17.3  kts.  Figure  9 also  shows  the  inertial 
tracks  of  the  ship  and  the  airplane.  The  aircraft  mean  bank  angle  therefore 
has  an  initial  bias,  maintained  during  the  transition  to  stationkeeping, 
since  direct  sideforce  control  is  not  available. 

6.2  Discussion  of  Results 

Figures  10(a),  (b)  and  (c)  show  the  principal  results  of  the  VOLAR  runs. 

The  dashed  lines  refer  to  limited  (saturating)  roll  and  yaw  controls.  The 
maximum  values  of  the  AV-8A  roll  and  yaw  control  moments  are  given  in  Section  3. 
The  solid  lines  represent  results  obtained  with  these  saturation  limits 
arbitrarily  removed,  i.e.  a completely  linear  ai rcraft  model . 

Figure  10(a)  shows  the  airwake  mean  and  random  components  as  predicted 
by  the  airwake  model  of  Reference  11.  The  magnitude  of  the  random  component 
follows  a generally  similar  trend  to  that  of  the  mean,  except  for  the  'spike' 

1 .i  ★ 

at  t = 7 secs.  The  physical  reason  for  this  spike  is  not  known,  but  it  is 
distinctly  snown  by  the  wind  tunnel  data  analysed  in  Reference  11  on  which 
the  airwake  model  of  Reference  12  is  based. 

Figure  10(b)  shows  the  mean  and  random  components  of  roll  and  yaw  control 
time  histories.  It  is  apparent  that  for  the  particular  flight  condition 
studied,  the  effect  of  control  limiting  is  small.  The  random  components  of 
the  control  time  histories  follow  a generally  similar  variation  to  the  random 
airwake  components.  Hence  there  is  considerably  more  random  control  activity 
at  the  beginning  of  the  approach  than  at  the  stationkeeping  stage,  where  the 
airplane  is  flying  above  the  boundary  of  the  airwake.  The  mean  roll  control 
deflection  displays  a spike  at  t = 7 secs  demonstrating  the  pilot's  reaction 
to  the  spike  in  the  mean  airwake. 

Figure  10(c)  shows  the  time  histories  of  bank  angle  and  lateral  position 
error.  The  latter  quantity  is  the  deviation  of  the  ground  track  from  the 
commanded  track.  The  mean  lateral  position  error  is  small  (<  1.2m)  but  the 
random  component  is  larger  (o  = 2.7m),  especially  at  the  beginning  of  the 


* 

On  Figure  10(a)  the  mean  and  random  components  are  designated  Vah  and  Vay 
respectively.  For  the  case  considered  here  these  terms  are  equivalent. 
The  commanded  flight  path  requires  the  airplane  to  fly  parallel  to  the 
W.O.D.  direction,  hence,  the  wind  side  component  referred  to  ship  wind 
axes  (Vay)  equals  aircraft  sideyust  component  (Vah) * 
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FIGURE  10(a).  AV-8A  APPROACH  WITH  CLASSICAL  PILOT  MODEL. 
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approach.  In  practice,  the  20  kt  ambient  wind/35  kt  W.C.D.  condition  would 
be  associated  with  appreciable  ship  motion  which  would  increase  the  error. 
Therefore  the  random  lateral  position  error  is  disturbing,  indicating  that 
the  AV-8A  lateral  characteristics  may  be  unsatisfactory.  Note  that  the 
bank  angle  time  histories  snow  a similar  ratio  of  mean  to  random  components, 
after  subtracting  the  initial  bias  required  to  follow  the  curved  ground 
track. 


7.0  OPTIMAL  PILOT  MODEL 


7.1  Background 

A disadvantage  of  tiie  classical  pilot  model  described  in  Section  5 is 
that  the  criteria  used  to  select  the  pilot  gains  are  based  on  verbal  adjust- 
ment rules  which  are  essentially  subjective.  To  remove  this  disadvantage 
Anderson  (Reference  ?l)*and  Dillow  (Reference  22)  replaced  tiie  verbal  adjustment 
rules  by  mathematical  optimization  criteria.  In  the  Anderson/Di  How  approach 
(known  as  “Paper  Pilot")  classical  loop  structures  are  retained  but  the  pilot 
gains  and  time  constants  are  found  as  the  result  of  minimizing  a performance 
index.  For  example,  for  the  longitudinal  hover  task  (with  no  coupling  into 
height  control)  the  performance  index  is  a linear  combination  of  r.m.s.  pitch 
rate,  r.m.s.  longitudinal  hover  error,  and  pilot  lead  time  constant. 

Paper  Pilot  type  models  have  been  developed  for  cruise  and  hover  control 
of  longitudinal  and  lateral  motions.  Each  of  these  four  cases  requires  a 
different  performance  index  and  a different  loop  structure.  It  was  originally 
proposed  to  use  the  Paper  Pilot  series  of  computer  programs  developed  by 
Anderson  and  Dillow  (References  22,  23,  2A)*as  subroutines  in  VOLAR.  Separate 
Paper  Pilot  programs  exist  for  longitudinal  hover,  lateral  hover,  longitudinal 
cruise  and  lateral  cruise.  A subsequent  publication  by  Dillow  (Reference  25)* 
replaces  these  separate  programs  by  a unified  model,  which  is  shown  to  give 
good  agreement  with  measured  data  from  simulator  tests,  (both  fixed-base  and 
moving-base).  Oil  low's  model  is  shown  in  Figure  11.  It  models  the  pilot  by  a 
Kalman  filter  in  series  with  an  optimal  controller.  The  optimal  controller 
minimizes  a performance  index,  J,  of  the  following  form: 

o - 

wnere 

y^  = the  ith  perceived  state  variable 

£ = rate  of  control  deflection 

Wy.R  are  weighting  functions 

Some  differences  between  this  index  and  the  performance  index  used  in  previous 
pilot  models  should  be  noted,  since  these  differences  appear  to  be  the  reason 


♦References  are  listed  on  page  50. 
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Figure  11.  Block  Diagram  of  the  Optimal  Pilot  Model 
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for  the  success  of  Dillow's  model.  For  example,  this  index  does  not  include 
all  state  variables,  only  those  perceived  by  the  pilot.  Note  also  that 
control  deflection  does  not  appear  in  the  index;  instead  control  rate  is 
employed.  As  explained  below  Dillow's  new  model  appears  to  be  better  suited 
to  our  purposes  than  are  the  original  paper  pilot  models. 

This  model  requires  some  generalization  to  the  mul ti -control  task  of 
V/STOL  Launch  and  Recovery  and  the  effect  of  airplane  flight  mode  (cruise, 
transition,  or  hover)  on  the  weighting  functions  demands  careful  study. 
Nevertheless  the  model  offers  the  following  advantages: 

(1)  A single  general  form  of  model  is  employed  for  all 
flight  conditions  (as  opposed  to  the  different  model 
forms  employed  for  "Cruise  Paper  Pilot".  Reference  23, 
and  "Hover  Paper  Pilot",  Reference  22). 

(2)  The  effects  of  visual  landing  aid  and  |_S0  cues 
can  be  incorporated  in  the  performance  index  by 
increasing  the  number  of  perceived  variables,  m,  in 
Equation  7.1  and  by  altering  the  weighting  functions. 

This  simplifies  modeling  the  LSO  and  visual 
landing  aids. 

(3)  One  general  formulation  is  employed  irrespective 
of  the  number  of  degrees  of  freedom  that  the 
pilot  is  controlling. 

On  the  other  hand,  some  disadvantages  of  the  optimal  controller  model  must 
be  noted. 

(1)  Relatively  little  multi-axis  data  has  been  analyzed 
to  yield  the  appropriate  performance  indices. 

(2)  An  estimate  must  be  made  (or  implied)  regarding  the 
relative  time  spent  controlling  each  axis,  for 
situations  where  the  axes  are  dynamically  uncoupled,  e.g., 
longitudinal  and  lateral  control.  Only  a small  data 

base  exists  for  such  estimates. 
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(3)  As  will  be  explained,  minimization  of  the  performance 
index  is  made  more  complicated  because  the  model 
assumes  that  the  covariance  matrices  describina 
measurement  and  neuromuscular  noise  are  functions  of 
the  numerical  value  of  the  performance  index. 

(4)  The  model  introduces  a large  number  of  additional 
states  because  of  the  estimator  and  predictor  dynamics. 

This  is  a significant  disadvantage  for  covariance 
propagation  programs  such  as  VOLAP  because  the  run 
time  of  such  programs  tends  to  vary  as  the  sauare  of 
the  number  of  state  variables.  (For  Monte  Carlo 
simulation  programs  run  time  tends  to  be  linearly 
proportional  to  the  number  of  state  variables.) 

7.2  The  Model 

This  section  reviews  the  evolution  of  the  optimal  pilot  model  from  the 
form  of  Figure  11  to  the  form  employed  in  VOLAP.  The  evolution  was  necessary 
to  balance  the  completeness,  accuracy,  and  efficiency  of  the  pilot  model 
versus  these  aspects  of  the  other  element  models  (shiD,  airwake,  etc.).  A 
general  discussion  is  presented  first,  followed  by  eguations  representing 
the  optimal  pilot  model  used  in  VOLAP. 

7.2.1  Simplification  of  Optimal  Pilot  Model 

The  model  presented  in  Figure  11  will  be  referred  to  here  as  the 
"complete  model".  This  form  of  the  optimal  pilot  model  has  been  validated  in 
References  25,  26,  27,  and  elsewhere. 
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27.  Stengel,  R.  F.,  J.  R.  Broussard,  P.  W.  Barry,  and  J.  H.  Taylor,  Modern 

Methods  of  Aircraft  Stability  and  Control  Analysis,  Report 
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Because  VOLAR  is  presently  programmed  only  for  continuous  systems,  the 
pure  time  delay  of  Figure  11  was  replaced  by  a first  order  Pade"  approxima- 
tion. This  removes  the  discrete  element  from  the  system  and  replaces  it  with  a 
continuous  one.  In  the  continuous  system  all  variables  are  defined  at  time  t. 
(For  the  system  of  Figure  11  to  be  optimal,  the  controller  must  operate  on 
x(t).  However  the  estimator  is  providing  estimates  that  are  shifted  back  in 
time  t seconds  because  of  the  pure  time  delay  on  the  measurements.  The 
predictor  is  required  to  furnish  (predict)  x(t)  to  the  controller  using 
x(t-T)  as  its  input.)  These  changes  are  reflected  in  the  model  shown  in 
Figure  12.  This  figure  shows  the  optimal  controller  pilot  model  that  we 
initially  employed  for  the  longitudinal  task.  At  the  time  of  its  use,  this 
pilot  model  was  part  of  a system  model  requiring  48  state  variables,  of  which 
30  were  attributed  to  the  pilot  model. 

Computer  run  time  is  approximately  proportional  to  the  square  of  the 
number  of  state  variables.  This  indicates  that  a simpler  pilot  model  would 
reduce  costs.  It  is  desirable  therefore  to  keep  the  number  of  state  variables 
low.  Besides  these  considerations  of  economy,  the  complete  pilot  model 
appears  unnecessarily  refined  for  the  present  study  when  compared  to  the 
state-of-the  art  ship  motion  and  airwake  models  employed. 

Eliminating  the  Pade  approximation  substantially  reduces  the  number 
of  state  variables  associated  with  the  pilot  model.  There  is  ample  justifi- 
cation in  the  literature  to  support  this  simplification.  For  instance. 
Reference  28  uses  the  optimal  pilot  without  the  Pade'  approximation  and  in 
Reference  29  the  authors  present  a strong  case  that  the  complete  model  is 
overparameterized,  i.e.  unnecessarily  complicated.  It  was  decided  to  drop 
the  Pad6  approximation  to  take  advantage  of  the  run  time  savings  and  bring 
the  pilot  model  more  in  line  with  the  size  and  complexity  of  the  other  models 
in  the  overall  system.  The  optimal  pilot  model  used  in  VOLAR  is  shown  in 
Figure  13  Note  tnat  in  Figure  13  the  block  depicting  the  neuromuscular  lag 


28.  Karmarkar,  J.  S.,  and  J.  A.  Sorensen,  Information  and  Display  Require- 

ments for  Independent  Landing  Monitors,  NASA  CR-2687,  1976. 

29.  Phatak,  A.,  II.  Weinert,  I.  Segall,  and  C.  N.  Day,  Identification  of  a 

Modified  Optimal  Control  Model  for  the  Human  Operator,  Automatica, 
Vol.  12,  1976. 


EXTERNAL  DISTURBANCES 


Figure  12.  Block  Diagram  of  Optimal  Pilot  Model  with  Pade  Approximati 


does  not  appear  explicitly.  An  effective  neuromuscular  lag  is  still  present 
but  is  implicit  in  the  optimal  controller.  This  is  achieved  by  altering 
the  weighting  functions  in  the  optimal  controller  performance  index,  as 
explained  in  Reference  26.  This  is  computationally  equivalent  to 
retaining  the  separate  block  representing  the  neuromuscular  dynamics.  No 
approximation  is  involved  bv  this  substitution. 

7.2.2  Equations  Describing  the  Optimal  Pilot  Model 

The  pilot  model  of  Figure  13  can  be  represented  by  linear  differential 
equations  with  time-varying  coefficients.  The  standard  Dilot  model,  (Eaua- 
tions  7.1  through  7.20)  without  the  time  delay  and  predictor  is  presented 
here  in  vector-matrix  format.  For  the  complete  set  of  eauations,  refer  to 
References  26  or  27. 

The  airframe  and  engine  dynamics  are  ni ven  by: 


x = Fj  x + Gj  u + Tw  1 
The  noise  coloring  processes  are  described  by: 


w'  = A w'  + B w 


(7.1) 


(7.2) 


where  w is  a vector  of  uncorrelated  Gaussian  white  noises  representing 
external  disturbances. 

The  controller  dynamics  for  the  airframe  and  engine  are  given  by: 


= C x + Rl  u + vu 


(7.3) 


where  vu  is  a vector  of  uncorrelated  Gaussian  white  noises  representing  pilot 
neuromuscular  ("motor")  noise.  The  elements  of  the  matrices  C and  R^  must  be 
determined  prior  to  using  VOLAR  and  are  the  solution  to  an  optimal  regulator 
problem.  The  determination  of  C and  R^  is  discussed  in  section  7.2.3  of 
this  report.  The  form  of  equation  7.3  is  standard  for  the  optimal  pilot 
model  formulation.  The  feedback  of  u into  u results  in  the  effective 
neuromuscular  lag  previously  mentioned.  A full  discussion  of  this  can  be 
found  in  Reference  26. 
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The  vector  y represents  the  measurements,  i.e.  the  quantities  perceived 

by  the  pilot.  The  first  derivatives  of  pilot  perceived  quantities  are  also 
included  in  y (References  25  and  26). 
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+ v. 


(7.4) 


The  quantity  v^  is  the  vector  of  measurement  noises.  These  are  Gaussian, 


uncorrelated,  and  white. 

The  estimator  dynamics*are  described  by: 
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where  the  matrix  K is  given  by: 
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The  value  of  PF  (the  estimator  covariance  matrix)  is  obtained  from: 


Pp  = 


<E  " PE 
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*A  good  general  discussion  of  estimators  is  given  in  Reference  30. 
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(7.12) 


The  covariance  matrix  for  the  estimated  states  is  obtained  from: 


P;  = P;  F*  + F-  P*  + PF  0 V"1  [0  0 0]  PF 

A A A A A L Jr  L 


(7.13) 


wnere: 


P*  = E j w'  w' 
/ u u 


(7.14) 
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(7.15) 


The  measurement  noise  covariance  matrix: 


[vy  vy] 


and  the  neuromuscular  noise  covariance  matrix: 


Vu  ■ E 


K vl] 


are  functions  of  tiie  covariance  matrix  Px: 
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where  Px  is  given  by: 


(7.16) 
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(7.17) 
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and  Vy  is: 


Vy  - Py  diag  ([D  0 0]  P, 
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(7.18) 


Noting  that  the  lower  right  hand  corner  of  Px  contains  a submatrix  pxu»w'iere: 


px  = E 
u 


{[»  uT)j 


(7.19) 


we  can  write: 

Vu 


= pu  diag  (Px  ) 


(7.20) 

Intuitively  one  would  expect 


The  cnoice  of  pu,  Py  requires  some  consent. 

Vu  and  Vy  to  be  constant.  However  Reference  26  suggests  that  Vu  and  Vy  are 


proportional  to  Px;  i.e 


pu  and  py  are  constants. 


More  experimental  data 
are  required  to  settle  this  point.  For  the  present  it  seems  advisable  to 
follow  the  approach  of  Reference  26,  for  which  supporting  data  exist. 

7.2.3  Calculation  of  C and  R. 


The  elements  of  C and  RL  in  equation  7.3  are  the  feedback  gains  of  the 
optimal  controller.  They  are  determined  by  solving  an  optimal  regulator 
problem.  In  the  equations  below  the  vector  u contains  aircraft  controllers 
(elevator,  throttle,  etc.)  and  uc  is  a vector  of  commanded  inputs. 


Given  the  system 
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determine  the  control  strategy  u*  that  minimizes  the  cost  function,  J. 


Note 


that  this  formulation  is  similar  to  the  classical  optimal  regulator  problem 
with  one  exception.  Here  the  conrol  vector  u has  been  adjoined  to  x and  w'. 
Observe  that  u = Iuc  in  equation  7.21  (I  is  the  identity  matrix);  therefore 
'control  rate'  is  commanded,  instead  of  'control  displacement'.  Consequently 
there  is  a weighting  matrix,  Rc>  on  control  rate  in  the  cost  function. 

The  solution  to  the  optimal  regulator  problem  is  given  in  Reference  30, 
as  follows: 

The  optimal  control  law  is 

(7-23) 
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Lhe  steady 
equation  given  below. 
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(7.24) 


and  P is  the  steady  state  solution  (i.e.  P = 0)  of  the  matrix  Riccati 
c t 
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Because  of  the  principle  of  separability  the  optimal  estimator  problem 
and  the  optimal  controller  problem  can  be  worked  independently.  Thus  one 
can  conveniently  solve  for  C and  R^  apart  from  VOLAR.  The  resulting  values 
of  C and  R^  become  inputs  to  VOLAR. 

7.2.4  The  Collected  Equations 

The  coupling  between  the  estimator  dynamics  and  the  system  dynamics  is 
described  by  equation  7.3 


u - C x + \ u + vu  (7.3) 


30.  Bryson,  A.  E.,  and  Y.  C.  Ho,  ApDlied  Optimal  Control,  Blairdell, 
Waltham,  Mass.,  1964. 
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This  equation  implies  that  the  estimator  equations  must  be  adjoined  to  the 
basic  system  equations  in  order  to  characterize  the  complete  pilot-plus- 
airframe  system.  Hence  if  the  system  dimensionality  is  of  order  n,  the 
pilot-plus-airframe  system  is  approximately  2n. 

The  overall  system  can  be  represented  as: 


(7.26) 


Figures  14  and  15  show  how  F and  G (respectively)  breakdown  into  their  sub- 
matrix components.  Note  that  for  simplicity  and  clarity  the  matrices  K and 
C have  also  been  partitioned. 

The  sub-matrix  breakdown  in  Figure  14  and  15  is  completely  general; 
i.e.,  it  is  valid  for  any  size  matrix.  The  actual  dimensions  on  each  of  the 
sub-matrices  are  problem-dependent.  Regardless  of  how  the  F matrix  is  parti- 
tioned it  must  contain  the  basic  system  dynamics  and  the  controller  and 
estimator  dynamics  representing  the  human  pilot. 

This  section  has  presented  the  equations  describing  the  present 
"optimal  pilot"  model.  Section  8 presents  an  example  demonstrating  the 
application  of  the  optimal  pilot  model.  Programing  details  for  using  the 
VOLAR  option  for  the  model  are  given  in  Appendix  A.  The  decision  whether  to 
use  the  optimal  pilot  model  or  the  classical  pilot  model  may  hinge  upon  the 
relative  run  times  required  by  the  VOLAR  program  for  each  mode.  Comparing  the 
"classical  pilot"  example  of  Section  6 with  the  "optimal  pilot"  example  of 
Section  8 yields  the  following  relative  run  times  for  a 20  second  time  history: 


Dimensionality  (order  of  F)  Run  Time  CPU  Secs  Relative  Run  Time 
Classical  17  54.4  1.0 

Optimal  39  355.9  6.5 

For  comparison,  a 200-run  Monte  Carlo  simulation  of  the  classical  pilot  example 
requires  800  CPU  seconds;  thus,  VOLAR  requires  only  7 percent  of  the  computer 
time  required  for  Monte  Carlo  analysis. 
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Figure  14  Elements  of  the  Generalized  F-Matrix 


Identity  Matrix 


Figure  15  Elements  of  the  Generalized  G-Matrix 


7.3  Modifications  for  Multi  Axis  Control 

Degradation  in  pilot  control  of  one  axis  due  to  an  increased  workload 
in  another  axis  effects  are  readily  included  in  the  optimal  pilot  model. 
These  effects  are  modeled  as  an  increase  in  the  noise  to  signal  ratio 
pyi. . Thus 


where 

py  = Noise  to  signal  ratio  associated  with 
"full  attention"  (Equation  7.18) 

fc  = Fraction  of  attention  devoted  to  the 
control  task  as  a whole 

f$  = Fraction  of  attention  devoted  to  sub- 
tasks; e.g.,  longitudinal  control 

f.j  = Fraction  of  attention  devoted  to  the 
ith  perceived  variable  and  its  first 
derivative  in  subtask  s;  e.g.,  control 
of  pitch  attitude  in  the  longitudinal 
subtask 

7.4  Visual  Landing  Aids 

Visual  landing  aids  and  Landing  Signal  Officer  (LSO)  functions 
can  be  easily  represented  by  appropriate  augmentation  of  the  pilot  perceived 
variables,  y,  in  Equation  7.4.  Pilot  perceived  variables  are  selected  to 
fit  the  task.  These  may  include  on  board  indicators  as  well  as  ship  based 
landing  aids.  Additional  perceived  state  variables,  such  as  the  LSO's  signals 
are  adjoined  to  y thereby  increasing  its  dimensionality.  If  the  LSO 
function  is  to  emphasize  certain  state  variables  which  are  already  included 
in  the  vector  of  pilot  perceived  variables  then  the  noises  on  those  variables 
must  be  adjusted  according  to  the  time  spent  viewing  each  of  the  sources. 

A discussion  on  the  treatment  of  scanning  behavior  can  be  found  in 
Reference  26. 
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8.0  AV-8A  LONGITUDINAL  APPROACH  IN  A LOW  SEA  STATE 

8.1  Task  Description 

The  simulated  task  is  a constant  altitude  decelerating  approach  to  a 
stationkeeping  position  50  ft.  (15.2  m)  above  the  deck,  followed  by  a 6 sec. 
descent  to  the  deck.  The  aircraft  is  an  AV-8A  flown  by  the  optimal  pilot 
model.  The  external  disturbances  modeled  are  Dryden  type  atmospheric  turbu- 
lence and  the  ship's  airwake.  Sea  State  2 is  assumed;  this  gives  only  small 
ship  motions,  which  are  considered  to  be  negligible  for  this  example.  (Ship 
motions  are  modeled  in  a subsequent  example  which  assumes  a higher  Sea  State.) 

Figures  16  and  17  illustrate  the  desired  approach  flight  path.  During 
the  initial  phase  the  airplane  decelerates  at  a constant  2.33  fps2  and  is 
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directed  to  the  ship  by  commanding  its  inertial  position.  The  ambient  wind 
speed  is  10.5  knots,  and  the  ship  is  steaming  at  15  knots.  As  shown  in 
Figure  16,  the  wind  and  ship  directions  are  such  that  the  W.O.D.  speed  is 
20  knots.  The  commanded  flight  path  requires  the  aircraft  to  be  headed  into 
the  W.O.D. , hence  the  commanded  sideslip  angle  is  zero.  It  is  therefore 
assumed  that  longitudinal  and  lateral  motions  are  decoupled,  and  only  longi- 
tudinal motions  are  simulated  in  this  example. 

8.2  Model  Components 

8.2.1  Airplane 

The  AV-8A  airplane  is  modeled  by  the  longitudinal  equations  of  motion 
discussed  in  Section  2.  The  dynamics  of  the  powerplant  are  modeled  as  a 
first  order  lag  with  a time  constant  of  0.2  seconds. 

The  airspeed  of  the  AV-8A  at  the  initiation  of  the  approach  segment 
simulated  here  is  24.6  knots  (17.81  m/sec)  and  the  terminal  airspeed  equals 
the  W.O.D.  speed.  AV-8A  flight  data  (Reference  35)  indicate  that  the  pilot 
flies  this  segment  of  the  approach  employing  only  two  longitudinal  controls, 
i.e.  thrust  magnitude  and  pitching  moment.  The  third  available  control, 
nozzle  angle,  is  set  at  a fixed  trim  setting.  Our  simulation  follows  the 
full-scale  example,  but  core  space  is  reserved  for  nozzle  angle  in  the  VOLAR 
program,  so  it  can  readily  be  included  as  a control  variable  if  desired. 

8.2.2  Atmosphere 

The  atmospheric  model  is  discussed  in  Section  4.  Dryden  turbulence 
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r.m.s.  values  for  ug  and  wg  were  1 ft/sec.  The  ship  airwake  was,  of  course, 
varied  as  a function  of  aircraft  position  with  respect  to  the  ship. 

8.2.3  Optimal  Pilot  Model 

A description  of  equations  representing  the  optimal  pilot  are  presented 
in  Section  7.  The  feedback  gains  of  the  optimal  controller  are  required  as 
input  to  VOLAR.  These  are  computed  separately.  This  is  a corresponding 
task  to  the  root  locus,  frequency,  and/or  time  history  analyses  that  are 
performed  'off  line1  when  using  the  classical  pilot  model.  The  work  of 
Dillow  and  Picha  (Reference  25)  was  selected  as  the  basis  for  the  optimal 
pilot  feedback  gain  calculation  used  in  VOLAR.  One  control  task  analyzed  in 
Reference  25  is  single  axis  longitudinal  position  control  for  a hovering 
aircraft.  This  task  is  somewhat  similar  to  the  terminal  phase  of  an  approach 
to  hover  followed  by  stationkeeping  with  the  ship.  A more  representative 
task  would  be  the  multi-axis  case  of  height  and  longitudinal  position  control. 
To  account  for  this  Dillow's  results  were  modified  slightly.  In  addition  to 
a pitch  controller,  a thrust  controller  was  added  and  the  weighting  on 
height  displacement  was  assumed  to  be  equal  to  the  weighting  on  longitudinal 
position  displacement.  Denoting  principal  diagonal  elements  of  Qc  in 
Equation  7.22  as  q-j-.  Dillow's  weightings  are  qg  = 0.4  sec2/deg2  and 
qv  = 1.0  ft"2.  The  added  weighting  for  height  control  is  qz  = 1.0  ft"2. 

The  elements  of  the  weighting  matrix,  Qc,  are  constant  for  a given  task, 
but  the  elements  of  the  control  weighting  matrix  Rc  vary.  Rc  is  a diagonal 
matrix  the  elements  of  which  are  varied  to  yield  a neuromuscular  lag  tine 
constant,  t^,  of  0.2  second  for  each  controller.  This  is  a standard  procedure 
wher,  using  the  optimal  pilot  model  as  described  in  Section  7.  For  this  AV-8A 
simulation,  at  20  knots  airspeed  the  values  of  the  main  diagonal  elements  of 
Rc  were  r^  = 0.0489  and  r6j  = 0.00603. 

The  numerical  solutions  for  the  feedback  gains  were  obtained  using  a 
computer  program  from  Reference  31.  This  program  was  selected  for  its 
availability;  however,  any  program  that  will  solve  the  matrix  Riccati  equation 
(Equation  7.25)  may  be  used. 


31.  Melsa,  James  L.,  and  Stephen  K.  Jones,  Computer  Programs  for  Computational 

Assistance  in  the  Study  of  Linear  Control  Theory,  McGraw-Hill,  Inc.,  1973. 
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Foil  owing  the  determination  of  the  optimal  controller  feedback  gains  the  user 
must  address  the  estimator.  Here,  as  before,  the  work  of  Dillow  and  Picha 
served  as  a guide.  The  values  of  py  and  pu  (Equations  7.18  and  7.20)  were 
taken  from  Dillow  and  are  O.OItt  and  0.003ir  respectively.  The  pilot  was 
assumed  to  perceive  the  variables  X)yP,  Z^yp,  and  e.  Standard  optimal  pilot 
model  practice  (References  25  and  26)  places  the  pilot  perceived  variables 
and  their  first  derivatives  in  the  measurement  vector  y (Equation  7.4).  This 
practice  was  followed  here  with  the  additional  assumption  that  the  pilot 
perceives  the  displays  (measurements)  simultaneously.  In  more  detailed 
analyses  (References  25  and  26)  the  pilot  is  assumed  to  divide  his  time 
between  the  various  displays.  This  results  in  some  degradation  in  pilot 
performance.  This  effect  could  be  modeled  by  VOLAR  if  desired,  but  would 
result  in  a pilot  model  which  is  over-sophisticated  relative  to  the  models 
currently  used  for  the  other  elements  of  the  Launch  and  Recovery  system. 

A similar  simplification  was  made  concerning  multi-axis  control.  The 
standard  optimal  pilot  model  apportions  the  time  spent  controlling  each  axis. 

In  other  words,  while  the  pilot  may  be  'flying'  all  axes,  he  mentally  gives 
more  priority  to  one  axis  over  another;  which  further  degrades  his  performance. 
Currently  the  optimal  pilot  in  VOLAR  does  not  include  this  additional  penalty. 
He  is  assumed  to  apply  his  full  attention  simultaneously  to  all  the  axes  that 
he  is  controlling.  This  assumption  is  justified  by:  (1)  the  degree  of  sophis- 
tication of  the  other  element  models  as  discussed  above  and,  (2)  the  lack  of 
experimental  data  on  multi -axis  task  control  apportionment.  Multi -axis 
control  apportionment  could  be  included  as  outlined  in  Section  7.3  as  more 
experimental  data  become  available. 


8.2.4  Ship  Motion 

The  sea  conditions  for  this  example  produce  negligible  ship  disturbances; 
i.e.  there  are  no  random  ship  motions.  The  only  ship  motion  is  the  mean 
motion  along  the  X earth  axis;  X^  = V<jHjp  cos  (V'shiP ^ * (T*16  Y earth  axis 

motion  is  not  relevant  in  this  longitudinal-only  simulation.) 


8.3  Simulation  of  Two-Stage  Approaches 

Current  practice  for  the  AV-8A  performing  an  approach  and  landing  to  a 
small  ship  involves  a two-stage  approach.  The  first  stage  consists  of  a 
decelerating  transition  from  some  initial  speed  and  altitude  to  a 'stationkeeping' 
flight  condition  where  the  aircraft  maintains  a constant  position  relative 
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to  the  ship.  (This  position  is  near  the  ship  but  may  be  outside  the  airwake.) 
The  second  stage  of  the  approach  comprises  the  final  descent  and  touchdown. 

This  stage  may  not  commence  until  an  appreciable  time  after  the  stationkeeping 
mode  has  been  established,  since  the  pilot  may  have  to  wait  for  deck  clearance, 
or  a lull  in  ship  motion.  The  duration  of  the  stationkeeping  interval  is 
therefore  highly  variable. 

For  the  present  simulation  the  stationkeepinq  interval  was  not  explicitly 
represented,  although  each  of  the  above  stages  was  simulated.  It  was  therefore 
necessary  to  do  two  runs,  one  for  each  stage,  and  to  re-initialize  the  state 
variables  at  the  commencement  of  the  second  stage.  This  re-initialization 
procedure  involves  re-setting  the  covariance  matrix P in  Equation  1.4.  It 
would  be  unrealistic  to  reset P equal  to  zero  because  at  the  commencement  of 
the  second  stage  the  aircraft  has  been  stationkeeping  for  an  unspecified 
period  in  the  presence  of  atmospheric  turbulence,  and  other  noise  sources. 

P was  therefore  obtained  by  solving  the  matrix  Riccati  equation  for  an 
abritrarily  chosen  time  t = 5 secs,  assuming  P = 0 at  t = 0.  No  problems 
were  encountered  with  this  technique  for  the  cases  studied  here  since  P 
converged  rapidly  to  a steady-state  value.  However,  for  such  initialization 
runs  it  is  desirable  to  employ  a higher  order  integration  algorithm  than  for 
the  rest  of  the  simulation.  Smaller  time  increments  should  also  be  employed. 
Detailed  comments  on  this  point  are  given  in  Appendix  A.  This  re-initializa- 
tion is  unnecessary  for  simulations  of  STOL  approaches,  VTOL  launches,  and 
other  continuously  defined  maneuvers,  including  VTOL  landings  that  are 
performed  without  a 'break'.  For  such  simulations  P is  initialized  only 
once,  to  provide  the  proper  initial  conditions  for  the  beginning  of  the 
interval  that  is  being  simulated. 

8.4  Results 

Figures  18(a)  through  18 ( i ) show  the  results  of  the  simulation  of  AV-8A 
longitudinal  approach  characteristics.  The  commanded  flight  profile  is  as 
shown  in  Figure  17  and  the  optimal  pilot  model  is  employed.  The  ship  motion 
is  negligible  for  the  assumed  Sea  State  2 condition.  The  figures  will  be 
discussed  below  in  three  groups,  as  follows: 

Figures  18(a)  and  (b)  - Airwake  Characteristics 

Figures  18(c)  and  (d)  - Control  Deflections  Generated  by  Optimal  Pilot 

Figures  18(e)  thru  (i)-  Aircraft  Response 
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Each  figure  shows  the  time  history  of  (i)  the  mean  and  (ii)  the  standard 
deviation  (a)  of  the  appropriate  quantity.  The  left  side  of  each  figure 
illustrates  the  decelerating  constant  altitude  stage  of  the  approach,  the 
right  side  illustrates  the  descent  stage  (the  left  and  right  sides  of  the 
figures  are  separated  to  indicate  that  the  second  stage  follows  the  first 
after  a stationkeeping  stage  of  indeterminate  duration).  For  the  present, 
only  the  dashed  lines  on  the  figures  should  be  studied.  The  solid  lines 
represent  a Sea  State  5 case  for  which  ship  motion  is  significant.  This  case 
is  discussed  in  Section  10,  and  compared  with  the  Sea  State  2 case. 

Airwake  Characteristics  - Figures  18(a)  and  (b)  show  the  airwake  horizontal 
and  vertical  velocity  components.  These  components  are  scaled  directly  from 
FF - 1052  wind  tunnel  data  as  described  in  Reference  11.  Note  that  the  mean 
and  standard  deviation  of  the  horizontal  component  becomes  most  intense  at 
the  beginning  and  end  of  the  approach.  A similar  trend  is  noted  for  the 
standard  deviation  of  the  vertical  component  aVz,  indicating  that  at  the 
stationkeeping  location  the  turbulence  level  is  low.  On  the  other  hand  a 
considerable  upwash  exists  at  this  location  (WAW  = -4.0  m/s). 

Control  Applied  by  Optimal  Pilot  - Figure  18(c)  shows  the  time  histories 
of  the  mean  and  standard  deviation  of  the  pitch  attitude  control,  06^.  The 
doublet  applied  by  optimal  pilot  to  mean  a6^at  t = 15  secs  is  required  to 
null  the  airplane  speed  relative  to  the  ship.  This  rapid  control  movement 
is  a mean  activity  only,  and  is  not  matched  by  a corresponding  increase  in 
random  stick  motions.  On  the  contrary,  the  latter  are  most  pronounced  at 
the  initiation  and  at  the  completion  of  the  recovery.  As  shown  in  Figure  18(d) 
the  throttle  time  histories  show  similar  characteristics  to  the  pitch  control 
motions . 

The  time  integrals  of  the  variances  of  the  random  control  deflections 
are  known  to  be  correlated  with  pilot  opinion,  as  shown  in  References  21,  22, 
and  25.  On  this  basis  one  would  expect  the  pilot  opinion  to  be  worse  at  the 
nr -ng  and  at  the  end  of  the  recovery,  when  the  airplane  is  immersed  in  the 
*'il*%*  turbulence  levels. 
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Aircraft  Response  - The  appropriate  pitch  attitude  time  histories 
are  shown  on  Figure  18(e).  As  one  would  expect,  a0  varies  with  time  in  a 
generally  similar  manner  to  the  standard  deviation  of  the  horizontal  airwake 
velocity  component.  The  maximum  magnitude  of  o0  is  approximately  0.6  degrees, 
indicating  that  for  this  sea  state  the  pilot  is  controlling  the  aircraft 
within  close  tolerances.  Figure  18(e)  indicates  relatively  large  variations 
in  mean  0 of  approximately  4°.  These  correspond  to  the  pull-up  maneuvers 
required  to  terminate  each  stage  of  the  approach.  These  large  terminal  pitch 
angles  could  be  trimmed  out  by  collective  nozzle  deflection,  but  as 
explained  previously  this  trimming  control  was  fixed  for  the  present  simulation. 

Figure  18(f)  shows  actual  and  commanded  altitude  time  histories.  The 
mean  height  error  is  very  small,  and  the  standard  deviation  is  less  than  0.2m 
for  this  case  with  negligible  ship  motion.  An  enlarged  plot  of  the  mean 
* height  error  is  given  on  Figure  18(g),  together  with  the  mean  longitudinal 

position  error.  The  latter  quanti ty  is  approximately  0.8m  at  touchdown,  with 
a corresponding  standard  deviation  of  approximately  0.25m  as  shown  in 
Figure  18(h).  As  shown  by  the  near-coincidence  of  the  lower  graphs  on 
Figure  18(h)  the  commanded  and  actual  positions  are  extremely  close  for  the 
final  part  of  the  approach. 

Correlation  of  the  above  results  with  flight  data  is  discussed  in 
Section  13.  On  the  basis  of  the  results  presented  here  the  recovery  at 
20  kts  W.O.D.  with  the  wind  direction  specified  in  Figure  16  appears  to  be  a 
maneuver  which  is  well  within  the  capabilities  of  the  pilot  and  the  aircraft. 
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9.0  SHIP  MOTION  MODEL 

9.1  Introduction 

The  examples  of  VSTOL  recoveries  presented  in  previous  sections  have 
been  for  low  sea  states,  where  ship  motion  is  nealiaible.  This  section 
describes  how  ship  motion  effects  are  modeled  by  VOLAP,  and  section  10 
presents  a numerical  example  for  a Sea  State  5 condition  for  which  shiD 
motion  is  important. 

Three  steps  are  required  to  include  ship  motion  in  the  state-variable 
equations  which  express  the  basic  covariance  propagation  relationship 
of  VOLAR. 

(1)  An  accurate  mathematical  model  must  be  found  for 
ship  response  to  specified  sea  states. 

(2)  The  model  must  be  approximated  bv  a model  which  is 
in  state-variable  form  and  which  does  not  unduly 
increase  the  dimensionality  of  the  state-variable 
equations  for  the  complete  svstem. 

(3)  The  approximate  model  must  be  adjoined  to  the 
equations  previously  formulated  for  the  case  of 
negligible  ship  motion. 

The  above  steps  are  described  in  Sections  9.2,  9.3,  and  9.4  respectively . 

9.2  The  RAOH  Ship  Motion  Model 

The  RAOH  ship  motion  model  is  described  in  Reference  32.  It  Dredicts 
ship  motion  spectral  characteristics  at  anv  aiven  freauencv  u rad/sec  by 
multiplying  wave  spectra  by  operators  known  as  "PA0"s  (response  amplitude 
operators).  For  the  example  of  Section  10  it  was  assumed  that  the  wave 
spectra  could  be  determined  from  the  Bretschneider  formula  (Reference  32), 


32.  Brown,  R.  G.,  and  F.  A.  Camaratta,  Navairengcen  Ship  Motion  Computer 
Program,  Naval  Air  Engineerina  Center,  Lakehurst,  New  Jersey, 
Report  NAEC-MISC-903-8,  1977.' 
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which  is  a function  of  the  significant  wave  height  (Hs)  and  the  modal  period 
(T0).  These  parameters  are  related  to  the  wind  and  shiD  velocity  vectors  as 
indicated  in  Table  8 (taken  from  Reference  33). 

For  the  present  application  the  RAO's  were  taken  from  data  supplied  bv 
NSRDC  for  a DD-963  ship. 

9.3  Approximation  of  the  RAOH  Ship  Motion  Spectra 

The  approximations  described  below  were  developed  by  P.  L.  Fortenbauoh 
in  Reference  33.  The  following  description  is  condensed  from  Reference  33, 
which  should  be  consulted  for  more  detailed  descriptions  of  modellinq 
procedures  for  ship-wave  interactions. 

The  numerical  data  for  the  ship  motion  spectra  calculated  by  RAO's  given 
in  Reference  32  were  approximated  by  equations  of  the  following  general  form. 

Gi  (jtoe)  = Ki  jaie 

f1'—)  * 2‘f  ^ (9.1) 

V *»n?  ) “"1 

where  i = the  motion  variable  = x,  y,  z,  es,  <}>s,  ip$ 

K.j  = Filter  gain 

uni  = The  equivalent  natural  freauenc.v  of  the  spectrum 
^ = The  equivalent  damping  ratio  of  the  spectrum 
u>e  is  the  wave  encounter  frequency 

The  form  of  Gi  (jwe)  was  determined  by  examination  of  the  spectra  developed 
for  the  conditions  in  Table  8.  For  the  most  part  these  spectra  assumed  the 
shape  shown  in  Figure  19  which  is  reasonably  represented  by  the  squared 
nagnitude  of  Gj  (j<ue).  It  is  expected  that  this  is  an  adequate  representation 


33.  Fortenbaugh,  R.  L.,  Application  of  the  NAFC  Ship  Motion  Simulation 
Program  for  Starboard  Approaches  to  00-963  Class  Ships,  Enclosure 
to  Vought  Report  2-55830/8AV0-153,  1978. 
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Table  8 Candidate  Ship  Speed-Wave  Direction-Wave 
Spectrum  Parameters-Wind  Over  Deck 
Combinations 


COMBIN- 

ATION 

SEA 

STATE 

vs 

y 

'Wd 

‘‘'WOD 

VWIND 

VW0D 

Hs 

To 

1 

6 

25 

120 

-60 

-30 

25.00 

43.30 

18 

15.13 

2 

5 

25 

120 

-60 

-30 

25.00 

43.30 

12 

13.50 

3 

5 

20 

120 

-60 

-30 

20.00 

34.64 

12 

13.50 

4 

5 

10 

135 

-45 

-30 

19.32 

27.32 

12 

13.07 

5 

5 

25 

180 

0 

0 

20*24 

45*49 

12 

12.07 

6 

5 

5 

180 

0 

0 

20*24 

25*29 

12 

11.51 

7 

4 

25 

105 

-75 

-30 

17.68 

34 . 15 

6.9 

10.6 

8 

3 

25 

105 

-75 

-30 

17.68 

34.15 

4.6 

8.8 

9 

3 

20 

105 

-75 

-30 

14.14 

27.32 

4.6 

8.8 

10 

3 

25 

90 

-90 

-30 

14.43 

28.87 

4.6 

8.8 

11 

3 

15 

120 

-60 

-30 

15.00 

25.98 

4.6 

8.8 

12 

3 

25 

180 

0 

0 

14*18 

39*43 

4.6 

8.8 

13 

3 

5 

180 

0 

0 

14*18 

19*23 

4.6 

8.8 

'‘'WIND  - 
^WOD  - 


Ship  Speed  (kt) 

Wave  Direction  (deg) 

Ambient  Wind  Direction  (deg) 
Wind  Over  Deck  Direction  (deg) 


WOD 


- Wind  Over  Deck  Soeed  (kt) 

- Significant  Wave  Weight  (ft) 

- Modal  Period  (sec) 
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for  all  single  mode  (having  one  well  defined  peak)  spectra  which  are  generally 
found  in  head,  bow,  and  beam  seas.  It  is  not  expected  to  apply  to  multimode 
spectra.  (One  of  the  combinations  - V$  = 5 kt,  p = 180°,  H$  = 4.6  ft, 

T0  = 8.8  sec  - had  a bimodal  spectrum  in  heave)  or  where  moderate  to  severe 
spectrum  warping  due  to  frequency  mapping  occurs  in  following  or  quartering 
seas. 

To  eastablish  Ki , , and  for  each  ship  degree  of  freedom  for  each 

condition  in  Table  8,  the  three  frequencies  Uemax»  WHI’  “LO^  ancl  maxi"mu,n 

gain  (K,  in  its  normal  units)  indicated  in  Figure  19  were  first  recorded. 

1 md  x 

The  two  frequencies  at  the  50  percent  peak  level  of  the  spectrum  establish 
the  damping  ratio  of  the  equivalent  spectrum  according  to  the  following 
relation: 


= (“HI  - “Lo) 
2uni 


(9.2) 


The  choice  of  the  50  percent  level  was  somewhat  arbitrary:  the  spectral  data 
indicated  that  generally  more  than  50  percent  of  the  total  power  was  contained 
in  the  - u>^q  frequency  band.  The  natural  frequency  and  damping  of  each 
filter  was  determined  as  follows.  It  was  assumed  that  and  wLq  are  the 
same  for  both  the  equivalent  and  measured  spectra.  The  natural  frequency 
for  this  case  is  given  by  “n, “ ^“HI  wL0  an<*  E9uation  9.2  again  used  to  calcu- 
late c.j.  Table  9 presents  the  results. 

Having  obtained  wn.  and  q was  determined  as  follows.  It  can  be 
shown  that  the  r.m.s.  output  of  G,-  (ju)0)  with  a white  roise  input  of  unit  r.m.s. 
is  given  by: 


Ki2  X3 


(9.3) 
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ENCOUNTER  FREQUENCY  (o>e)  * RAD/SEC 
Figure  19.  Typical  Ship  Motion  Response  Spectrum 
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Using  oi  for  the  32  frequency  component  approximation  data  of  Reference  32, 
Equation  (4)  was  solved  to  give  : 


The  resulting  K^'s  are  given  in  Table  9. 

Application  of  the  equivalent  filters  defined  in  Table  9 to  represent 
ship  motion  spectra  produces  only  the  autocorrelation  (power)  spectra.  The 
cross  spectra  require  further  study  and  definition.  In  terms  of  the 
covariance  matrix  for  ship  motions,  the  main  diagonal  contributions 
(variances)  have  been  defined,  the  off-diagonal  contributions  (covariances) 
have  not.  Knowledge  of  these  cross-correlation  effects  would  increase  the 
accuracy  of  the  simulation.  However,  in  view  of  the  numerous  approximations 
that  have  been  made  to  obtain  the  power  spectra  data,  it  appears  desirable 
to  defer  any  attempt  to  define  cross-spectra  until  a coordinated  improvement 
of  the  accuracy  of  all  system  element  models  is  made.  At  such  time  the 
ship  motion  power  spectra  models  should  also  be  refined  and  updated. 

9.4  Incorporation  of  Ship  Motion  Model  in  VOLAR 

Let  us  consider  the  simulation  of  the  longitudinal  ship  motions,  surge, 
heave,  and  pitch.  The  spectra  given  in  the  previous  section  all  conform  to 
Equation  9.1.  Because  of  this,  and  also  because  the  ship  surge,  heave,  and 
pitch  motions  are  correlated  in  an  unspecified  manner  they  may  be  simulated 
by  three  parallel  filters  forced  by  a common  white  noise  generator  as  shown 
in  Figure  20.  (If  the  motions  had  zero  correlation,  three  separate 
uncorrelated  white  noise  sources  would  be  required.  If  the  cross  spectra 
had  been  specified,  additional  separate  white  noise  generators  would  be 
required  with  appropriate  cross-filter  circuits.) 

In  Figure  20,  the  variables  e^,  e2»...eg  denote  additional  state 
variables.  The  appropriate  state-variable  equation  for  Figure  20  is  given 
on  page  88.  This  equation  must  be  adjoined  to  the  equations  given  previously 
for  the  case  of  zero  ship  motion. 
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NOTE:  *BI MODAL  RESPONSE  - DATA  CORRESPOND  TO  PEAK  WIDTH  SELECTED  TO  FIT  TOTAL  PESPONSE  WITH  ONE  PEAK 


Equation  for  Simulation  of  Correlated  3 DOF  Ship  Motion 


10.0  EXAMPLE  OF  AV-8A  RECOVERY  IN  A HIGH  SEA-STATE 


10.1  Task  Description 

The  simulation  presented  in  this  section  represents  an  AV-8A  performing 
an  approach  maneuver  which  is  similar  to  that  given  in  Section  8 except  that 
the  ship  motion  is  not  negligible.  A higher  Sea  State  is  selected;  Sea 
State  5 is  now  assumed.  Ship  motion  is  included  in  the  VOLAR  simulation 
using  the  formulations  described  in  the  preceding  section. 

Some  changes  must  be  made  in  the  ship  inertial  speed  and  ambient  wind 
to  assure  compatibility  of  the  airwake  with  the  increased  sea  state.  As 
illustrated  in  Figure  21  the  ship  inertial  speed  is  20  knots  at  a heading 
of  30  degrees.  The  ambient  wind  speed  is  20  knots  at  a heading  of  150 
degrees,  yielding  a W.O.D.  of  35  knots  at  0 degrees  heading  relative  to  the 
commanded  aircraft  heading.  Thus  the  aircraft  is  flying  into  the  W.O.D.,  as 
in  the  example  given  in  Section  8.  However,  in  that  example  a Sea  State  2 
and  a lower  W.O.D.  were  assumed. 

10.2  Element  Models 

The  AV-8A  airframe,  control  system,  and  powerplant  models  are  the  same 
as  those  used  in  previous  sections,  but  are  extended  to  higher  airspeeds. 

The  rate-of-closure  with  the  ship  is  the  same  as  in  the  example  of 
Section  8,  but  the  inertial  speeds  and  airspeeds  have  been  changed  because 
of  differences  in  ship  speed  and  ambient  wind  speed.  Keeping  rate-of- 
closure  constant  in  the  two  examples  seems  logical.  A pilot  in  a VTOL 
approach  may  not  be  aware  of  wind  speed  or  ground  speed  but  he  can  perceive 
rate-of-closure  with  the  ship.  He  can  fly  the  airplane  with  this  parameter 
fixed  and  accept  whatever  airspeed  results  since  a VTOL  aircraft  is  not 
limited  to  some  minimum  airspeed  by  stall  characteristics. 

For  this  simulation,  the  initial  airspeed  is  44.7  knots  (23  m/s)  and 
the  terminal  airspeed  equals  the  W.O.D.  speed.  As  in  the  previous  simulation, 
the  pilot  employs  two  longitudinal  controls;  thrust  magnitude  and  pitching 
moment.  It  is  assumed  that  the  nozzle  angle  is  not  re-triitmed  during  the 
approach.  The  Dryden  turbulence  spectra  r.m.s.  values  are  the  same  as  in 
the  previous  example  ug  = wg  = 1 ft/sec.  The  ship  airwake  is  different 
because  of  the  higher  W.O.D. 
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Figure  21  Commanded  Approach  Profile  for  Sea-State  5 Condition 
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The  optimal  pilot  feedback  gains  were  determined  as  discussed  in 
Section  8 for  two  points  along  the  trajectory.  The  values  for  intermediate 
airspeeds  were  obtained  by  linear  interpolation.  Note  that  the  elements  of 
the  weighting  matrix  Qc  (Equation  7.22)  are  the  same  here  as  in  Section  8 
because  they  are  task-dependent.  The  feedback  gains  are  different  because 
the  dimensional  derivatives  vary  with  airspeed.  The  remaining  details  of 
the  pilot  model  are  as  described  in  Section  8. 

10.3  Ship  Motion 

Two  problems  arise  in  including  the  ship  motion  in  a simulation  of  the 
longitudinal  degrees  of  freedom  of  a launch  or  recovery.  Firstly,  the  air- 
plane plane  of  symmetry  is  not  parallel  to  the  plane  of  symmetry  of  the  ship. 
Therefore  ship  motion  components  (heave,  surge,  sway,  pitch,  roll,  yaw)  that 
are  conventionally  referred  to  ship  centerline  axes  must  be  transformed  to 
an  axis  system  parallel  to  the  airplane  body  axis  system.  The  second  problem 
is  that  it  is  not  clear  to  what  extent  the  pilot  tracks  the  ship  motion. 

In  the  first  stage  of  a visual  approach  he  may  decide  to  ignore  ship  motion 
since  it  does  not  become  relevant  to  his  safety  until  after  he  terminates  the 
stationkeeping  mode  and  descends  to  the  deck.  The  ability  of  the  pilot  to 
detect  ship  motion  is  also  not  known.  Thresholds  must  exist  for  pilot 
detection  of  each  ship  motion  degree  of  freedom,  but  these  thresholds  may 
depend  on  visual  cues  from  superstructure  geometry  which  may  vary  from  one 
class  of  ship  to  another. 

To  demonstrate  the  VOLAR  program  it  was  assumed  that  the  pilot  tracks 
the  ship  heave  motion,  smoothed  by  a first-order  filter  as  shown  in  Figure  22 
This  implies  that  the  pilot  is  primarily  concerned  about  the  low  frequency 
components  of  the  ship  heave  motion.  The  only  ship  degree  of  freedom 
simulated  is  heave.  Some  of  the  other  ship  degrees  of  freedom  are  known  to 
be  small  (e.g.  surge),  and  the  cross-spectral  (covariance)  data  required  to 
properly  simulate  the  phasing  between  ship  motions  were  not  available.  For 
example,  computation  of  the  heave  motion  of  a general  point  on  the  ship's 
axis  of  symmetry  requires  the  heave  spectrum,  the  pitch  spectrum,  and  the 
pitch-heave  cross-spectrum.  When  such  cross-spectral  data  become  available 
they  can  readily  be  Incorporated  within  the  existing  VOLAR  input.  Without 
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Figure  22  Modified  Pilot  Altitude  Command 
due  to  Ship  Heave  Motion 
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such  cross-spectral  data  there  is  no  point  in  performing  the  ship/airp'.ane 
axis  transformation  mentioned  above.  By  considering  only  heave  ship  motion 
this  problem  is  avoided  in  the  present  example. 

10.4  Results  of  the  Example 

The  VOLAR  results  are  given  on  Figure  18,  (page  69),  together  with  results 
for  the  case  with  negligible  ship  motion  (Sea  State  2).  The  latter  results  have 
already  been  discussed  in  Section  8.  As  before,  the  stationkeeping  stage 
is  not  shown;  only  approach  to  stationkeeping  and  descent  are  presented 
in  Figure  18- 

Airwake  (Figures  18c,  d)  - Referring  first  to  Figure  18(b)  it  is 
apparent  that  the  increase  in  W.O.D.  speed  from  20  knots  to  35  knots  increases 
the  r.m.s.  turbulence  levels  in  the  airwake.  Although  a similar  general 
trend  exists  for  the  mean  components  it  is  interesting  to  note  that  for 
7 < t < 20  secs  the  u-component  of  the  mean  airwake  is  smaller  for  the 
35  knot  W.O.D.  case  than  for  the  20  knot  condition.  This  phenomenon  exists 
in  the  basic  wind  tunnel  data  analysed  in  Reference  11.  No  explanation  is 
apparent.  In  view  of  the  reasonable-looking  nature  of  the  remainder  of  the 
airwake  data  it  was  decided  not  to  adjust  this  apparently  anomalous  character- 
istic to  fit  preconceived  ideas  of  airwake  characteristics. 

Control  Deflections  (Figures  18c,  d)  - These  figures  indicate  a large 
increase  in  random  control  activity  between  the  20  knot  and  35  knot  W.O.D. 
cases.  The  one-sigma  random  throttle  movements  increase  approximately  by  a 
factor  of  6.  The  corresponding  factor  for  the  pitch  control  movements  is 
approximately  3.  The  difference  in  these  factors  is  expected.  The  difficulty 
of  the  pilot's  height  control  task  has  increased  not  only  because  of  increased 
turbulence  but  also  because  of  the  ship  heave  motion,  whereas  his  longitudinal 
control  task  has  become  more  difficult  only  because  of  increased  turbulence. 

Aircraft  Response  - Figure  18(e)  shows  the  p'tch  attitude  time  history. 

The  mean  e is  increased  compared  to  the  20  knot  condition.  The  high  mean 
pitch  altitudes  shown  result  partly  from  the  fact  that  the  thrust  nozzle 
vector  control  was  not  used. 

Figure  18(f)  illustrates  the  random  component  of  the  height  response. 

There  is  a marked  increase  in  this  quantity  associated  with  the  higher  sea 
state.  This  increase  occurs  despite  the  increased  thrust  control  activity 
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shown  or,  Figure  18(d),  indicating  that  the  pilot  is  more  active  but  less 
successful  in  controlling  the  inertial  altitude  of  the  airplane. 

As  the  instant  of  touchdown  is  approached,  inertial  altitude  becomes 

less  significant  than  the  altitude  of  the  airplane  relative  to  the  heaving 

deck.  The  mean  and  random  values  of  this  quantity  are  shown  in  Figures 

18(g)  and  (h)  respectively.  Figure  18(g)  illustrates  the  mean  vertical 
★ 

position  error,  which  is  less  than  0.5m  for  both  the  sea  states  simulated. 
Figure  18(h)  shows  the  random  component  of  the  airplane  altitude  relative  to 
the  ship.  The  lo  value  of  this  quantity  is  approximately  0.6m  throughout 
the  recovery.  This  is  a marked  increase  above  the  Sea  State  2 case  where  the 
corresponding  value  was  less  than  0.1m. 

The  longitudinal  position  error  is  shown  on  Figure  18(g).  Some  degrada- 
tion in  performance  is  evident  relative  to  the  Sea  State  2 case;  the  mean 
position  error  at  touchdown  increases  from  0.5m  to  1.5m.  The  corresponding 
random  component,  shown  on  Figure  1 8 ( i ) , shows  a similar  degradation  during 
the  first  stage  of  the  approach  but  not  at  touchdown.  This  is  apparently 
due  to  the  anomalous  airwake  model  characteristics  discussed  at  the  beginning 
of  this  section. 

The  lower  part  of  Figure  18(i)  shows  commanded  and  actual  longitudinal 
position.  This  shows  that  the  aircraft  acquires  the  commanded  approach 
profile  without  any  oscillations. 

In  summary,  the  example  results  presented  here  indicate  that  the  predicted 
behavior  of  the  aircraft-ship-pilot  system  is  consistent  with  the  assumptions 
made  for  modeling  each  element.  A comparison  with  flight  tests  is  given  in 
Section  13. 


*The  mean  vertical  position  error  shown  on  Figure  18(g)  is  the  mean 
difference  between  the  airplane  inertial  altitude  and  the  commanded 
altitude,  which  as  explained  in  Section  10.3,  includes  a filtered 
component  of  ship  motion.  This  filter  is  not  included  in  the  quantity 
o(ZAP  - Zs)  shown  in  Figure  18(h). 
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11.0  PROXIMITY  EFFECTS 
11.1  List  of  Proximity  Effects 

Table  10  lists  several  distinct  phenomena  which  cause  the  forces  and 
moments  acting  on  an  airplane  to  depend  on  the  proximity  of  the  airplane  to 
the  ship.  Only  very  limited  data  are  available  on  most  of  these  proximity 
effects,  as  indicated  below.  (The  numbers  refer  to  Table  10.) 

1. A  Reference  14  presents  extensive  look-up  tables  for  STOL-type  power- 

off  ground  effect.  These  tables  were  apparently  computed  on  the 
basis  of  an  Infinite  ground  plane.  As  shown  in  Reference  34  this 
assumption  may  be  invalid  for  aircraft  operating  from  ships, 
particularly  so  for  small  ships  where  the  ratio  of  span  to  deck 
width  Is  large  (i.e.  > 0.5). 

2.  A Fountain  and  suckdown  effects  are  highly  configuration-dependent. 

Reference  13  models  the  AV-8A  dominant  suckdown  effect  as  a 
6 percent  thrust  loss  at  touchdown,  decreasing  to  zero  at  a deck 
clearance  of  7 meters.  Reference  14  presents  a fairly  complex 
model  for  the  AV-8B  fountain/suckdown  effect.  The  AV-8A  and 
AV-8B  have  different  fountain/suckdown  characteristics  because  the 
AV-8B  has  a fountain  box. 

2.B  Reingestion  is  the  process  whereby  exhaust  gases  diffuse  and  then 
re-circulate  around  the  aircraft  due  to  convection  and  relative 
wind.  The  AV-8A  NATOPS  flight  manual  (Reference  35)  indicates 
that  reingestion  is  a potentially  serious  problem  for  protracted 
VSTOL  operations  in  which  the  airplane  remains  within  the  same  air 
mass  for  sufficient  time  for  re-circulation  to  develop.  Reference 
35  also  notes  that  a more  rapidly  acting  type  of  reingestion 
occurs  during  STOL  operations  due  to  part  of  the  forward  jet 
exhausts  impinging  on  the  ground  plane  and  being  deflected  forward 
just  ahead  of  the  Intake  where  the  gases  are  reingested.  Apart 
from  qualitative  comments  (e.g.  Reference  35)  no  AV-8A  reingestion 
data  have  been  published. 

34.  Levey,  H.  C.,  The  "Ground"  Interference  of  a Carrier  Deck, 

J.  Roy.  Aero.  Soc.,  April  1957,  pp.  276-281. 

35.  NATOPS  Flight  Manual,  NAVAIR  Rept  01-AV8A-1,  1976. 
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TABLE  10  PROXIMITY  EFFECTS 


1.  "Ground  Effects" 

2.  "Ground  Effects" 

mrnrnm 

4.  Restraining 

on  aerodynamic 

associated  with 

reactions 

system  forces 

forces  and  moments 
of  the  power-off 
configuration 

1A.  Ground  plane 

liftinq  enqine 
thrust 

2A.  Fountain  effects 

3A.  These  include 

4A.  Haul -down 

alters  induced 

and  suckdown  alter 

brake  forces, 

and  securing 

angle  of  attack 

mean  forces  acting 

landing  gear 

gear. 

of  lifting 

on  the  airplane. 

reactions,  and 

surfaces . 

oleo  damping. 

4B.  Launch  and 

2B.  Reingestion 

The  latter 

arrester  gear 

effects  are  slowly 

modify  the  ship 

where 

time-varying. 

motions  trans- 

appropriate. 

2C.  The  random  forces 
are  also  affected 
since  the  aircraft 
is  immersed  in 
self-generated 
turbulence. 

2D.  Changes  in  these 
effects  are  induced 
by  irregular  ground 
planes  (e.g.,  "edge 
of  deck",  or 
"sculptured  deck"). 

mitted  to  the 
airplane. 

3B.  Ski -jump 
effects . 

i 

2C.  AV-8A  flight  experience  (Reference  35)  indicates  that  the  aircraft 
self-generated  turbulence  is  a significant  factor  in  ground  effect. 

No  quantitative  data  have  been  published  on  this  phenomenon  for 
the  AV-8A. 

2D.  Edge-of-deck  effects  are  expected  to  be  important  for  VTOL  configura- 
tions in  conditions  where  fountain  effects  are  large  (e.g.,  AV-8B 
rather  than  AV-8A).  No  quantitative  data  have  been  published  for 
VSTOL  configurations  other  than  helicopters,  although  model  tests 
are  currently  being  planned. 

3A.  Data  on  AV-8A  and  AV-8B  landing  gear  reactions  are  given  in 
References  13  and  14. 
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3B.  The  "ski-jump"  STO  technique  (Reference  34)  involves  a curved 
deck  which  acts  through  the  landing  gear  reactions  to  produce  an 
upward  acceleration  for  STO.  If  data  are  supplied  regarding  the 
landing  gear  reactions  and  the  deck  shape,  the  applied  acceleration 
can  be  computed  without  difficulty. 

4A.  Haul-down  and  securing  systems  for  helicopters  are  described  in 
Reference  36.  No  data  have  been  published  on  such  systems  for 
the  AV-8A  or  AV-8B. 

11.2  Simulation  of  Proximity  Effects 

The  above  summary  shows  that  there  is  inadequate  data  on  many  proximity 
effects  to  permit  accurate  models  to  be  constructed.  This  is  particularly 
true  for  random  or  unsteady  components  such  as  self-generated  turbulence. 

For  covariance  propagation  simulations  the  mean  and  random  components  may  be 
equally  significant  and  it  is  not  justifiable  to  simulate  one  without  the 
other  without  experimental  evidence  to  justify  such  an  approximation. 
Accordingly,  proximity  effects  were  not  included  in  the  examples  presented 
In  this  report.  When  adequate  data  become  available,  proximity  effects  may 
be  modeled  as  indicated  in  Figure  23.  This  figure  illustrates  the  data 
processing  required  for  simulation  of  the  aerodynamic  proximity  effects 
listed  in  columns  1 and  2 of  Table  10.  Simulation  of  landing  gear  reactions 
may  be  accomplished  as  in  Reference  13,  and  restraining  system  forces  can  be 
included  with  the  landing  gear  reactions  if  desired. 


36.  Fozard,  J.  W. , Sea  Harrier  - The  First  of  the  New  Wave,  Aeronautical  J., 
January  1977,  pp.  15-40. 


37.  Bryson,  L.  B.,  F.  E.  Heenan,  C.  A.  Johnson,  Helicopters  In  the  Royal 
Navy,  Aeronautical  J.,  August  1972,  pp.  469-500. 
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FIGURE  23  FLOW  CHARI  FOR  SIMULATION  OF  PROXIMITY  EFFECTS 
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12.0  MISCELLANEOUS  TOPICS 


12.1  Introduction 

This  section  discusses  some  topics  which  do  not  fit  under  any  of  the 
previous  section  headings.  These  include  the  use  of  probability  ellipses 
to  present  the  mean  and  variance  Information  on  a single  graph  and  the 
simulation  of  launch,  in  addition  to  recovery. 

12.2  Probability  Ellipses 

The  results  of  the  examples  presented  earlier  in  this  report  have  been 
presented  as  graphs  of  the  mean  and  standard  deviation  of  appropriate  state 
variables  versus  time.  A disadvantage  of  this  type  of  presentation  is  that 
two  graphs  are  required  to  Illustrate  the  time  history  of  each  state  variable, 
one  for  the  mean  and  one  for  the  standard  deviation  of  the  random  component. 
This  disadvantage  can  be  overcome  by  using  the  form  of  presentation  sketched 
in  Figure  24.  This  consists  of  a graph  showing  the  mean  values  of  two  state 

variables  (or  output  variables)  Xj  and  x£»  with  the  mean  value  denoted  by 

^2"  As  shown  Reference  30,  the  covariance  matrix  P.of  Xj  and  ^ can 

be  related  to  the  probability,  p,  that  the  vector  x with  components  Xj 

lies  within  an  ellipse  defined  by  the  equation 

(x  - x)T  P_1  (X  - Y)  = L2  (12.1) 

For  jointly  Gaussian  distributions  of  Xj  and  x2  the  above  probability,  P,  is 
given  by 


p ■ 1 - e ('L2/2ax1  °x2)  (12.2) 

(fcxj  °x2) 

The  ellipse  axes  are  tilted  because  (in  general)  P is  not  a diagonal  matrix. 
In  other  words  the  covariance  or  cross-correlation  of  xi  and  x2  causes  the 
ellipse  axes  to  be  nonparallel  to  the  Xj  and  x2  axes. 
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FIGURE  24  COMPACT  PRESENTATION  SHOWNING  MEAN 
VALUES  AND  DEGREE  OF  RANDOMNESS  OF 
TWO  VARIABLES. 


In  principle,  the  form  of  presentation  shown  In  Figure  24  Is  more 
compact  than  plotting  xj,  X2»  oXj,  o versus  time  on  four  separate  graphs, 
as  in  the  examples  presented  earlier.  VOLAR  includes  the  capability  to 
output  probability  ellipses.  Figure  25  shows  a typical  example  generated 
by  VOLAR.  This  example  illustrates  a practical  problem  that  can  occur  with 
this  type  of  presentation.  The  problem  relates  to  scaling.  It  Is  often 
difficult  to  pick  comnon  scales  for  xj  and  oXj,  and  for  ^ and  that  will 
not  result  In  the  ellipses  becoming  too  small  or  too  narrow.  One  can  of 
course  choose  L to  yield  a large  probability  thus  enlarging  the  ellipse,  or 
one  can  employ  a different  scale  for  oXl*  retaining  the  original  scales  for 
7j,  X2»  and  ox The  latter  alternative  changes  the  aspect  ratio  of  the 
ellipse.  It  is  a matter  of  user's  choice  whether  this  distortion  is 
preferable  to  the  lack  of  resolution  implied  by  a very  narrow  but  unsealed 
ellipse.  In  Figure  25  both  the  o-axes  are  drawn  to  5 times  the  common  scale 
of  the  x-axes,  distorting  the  scale  of  the  ellipse  but  not  Its  proportions. 

12.3  Simulation  of  Launch 

In  the  examples  presented  in  this  report  VOLAR  has  been  employed  to 
simulate  recoveries  In  various  Sea  States,  but  no  launches  have  been  simulated. 
The  emphasis  on  recovery  reflects  flight  experience  that  this  Is  a more 
demanding  task  than  launch.  VOLAR  may  be  applied  to  launch  and  Indeed  to 
any  other  maneuver  that  can  be  simulated  by  conventional  Monte  Carlo  methods. 

To  simulate  launch  It  is  suggested  that  the  landing  gear  dynamics  should 
be  neglected,  and  the  aircraft  position  and  velocity  merely  initialized  at 
deck  levels.  If  a more  precise  simulation  is  desired  it  will  be  necessary 
to  model  the  landing  gear  dynamics. 

Currently  VOLAR  does  not  Include  a subroutine  for  landing  gear  reactions. 
Such  a subroutine  could  be  added  using  the  Monte  Carlo  program  of  Reference  13 
as  a guide.  It  would  be  desirable  to  formulate  the  landing  gear  reactions  In 
a form  that  Is  sufficiently  general  to  model  ski-jump  type  launches.  The 
comments  of  Section  11  regarding  current  deficiencies  In  analytic  models  of 
proximity  (ground)  effects  should  also  be  noted  since  such  effects  may  be 
Important  during  launch  operations. 
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Figure  25  Example  VOLAR  Printout  of  Probability 
Ellipses 
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13.0  CORRELATION  OF  SIMULATED  APPROACHES  WITH 
FLIGHT  TEST  RESULTS 


13.1  Introduction 

In  this  section  the  results  of  the  VOLAR  simulations  of  Section  6,  8, 
and  10  are  compared  with  flight  data  on  the  AV-8A.  A quantitative  data  base 
for  the  flight  tests  is  given  in  Reference  9 which  reports  AV-8A  shipboard 
suitability  trials.  In  these  trials  the  operating  envelope  of  the  AV-8A 
was  determined  partly  by  performing  VTO's  and  VL's  but  primarily  by  attempting 
to  hover  (i.e.  stationkeep)  over  all  portions  of  the  deck  of  an  LPD-class 
ship  while  varying  the  W.O.D.  as  much  as  possible.  A boundary  was  established 
relating  combinations  of  the  W.O.D.  magnitude  and  W.O.D.  angle  for  which  the 
pilot  was  able  to  stationkeep.  This  boundary  is  plotted  on  a polar  graph 
with  W.O.D.  magnitude  as  the  radial  coordinate  and  W.O.D.  direction  relative 
to  the  ship  centerline  as  the  angular  coordinate.  Figure  26  shows  this 
graph,  which  is  known  as  a "dynamic  interface  chart".  Figure  26  is  taken 
from  Reference  9. 

In  Figure  26  the  shaded  area  represents  the  combination  of  W.O.D.  magni- 
tude and  direction  for  which  the  above  described  AV-8A  operations  were  found 
to  be  feas  ble  for  starboard  approaches.  The  deck  outline  is  indicated  by  the 
heavy  lines,  and  the  deck  markings  for  the  two  landing  pads  are  represented  by 
lines  of  lighter  weight.  The  direction  of  the  approach  flight  path  relative 
to  the  ship  is  not  necessarily  parallel  to  any  of  these  deck  lines.  These 
lines  are  a useful  reference  to  the  pilot  and  may  be  used  by  him  to  establish 
a safe  approach  direction,  but  they  do  not  indicate  an  "optimum"  approach 
direction.  (Even  if  such  an  optimum  were  known,  it  would  vary  with  W.O.D. 
direction  and  hence  could  not  be  consistently  indicated  by  fixed  lines  painted 
on  the  deck.)  This  point  is  emphasized  because,  for  the  examples  in  this 
report,  the  approach  was  made  at  150  degrees  to  the  ship  centerline,  with 
the  flight  path  parallel  to  the  wind  over  deck.  As  shown  in  Figure  26  this 
approach  is  approximately  15  degrees  to  port  of  the  appropriate  deck  line. 

The  heavy  dots  in  Figure  26  denote  the  two  flight  conditions  simulated  by 
VOLAR.  The  20  knot  W.O.D.  condition  discussed  in  Section  8 is  within  the 
capabilities  of  the  pilot-airframe-ship  system.  The  35  knot  conditions 
discussed  in  Sections  6 and  10  lies  outside  the  operational  boundary. 
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EXAMPLE  CASE  OF 
SECTION  10 


EXAMPLE 
CASE  OF 
SECTION  8. 
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Figure  26  LPD  Dynamic  Interface  Chart 


The  type  of  AV-8A  used  in  the  flight  tests  of  Reference  9 was  an  early 
version  of  the  MK.  50  Harrier  - one  of  the  first  twelve  AV-8A'$.  These 
aircraft  had  only  a roll  stability  augmenter  system  which  was  of  very 
limited  authority,  no  head-up-display  was  employed,  and  only  daylight 
landings  were  performed.  (Night  landing  attempts  were  made  but  these  were 
not  successful.)  The  AV-8A  model  employed  for  the  VOLAR  simulations  did 
not  include  any  S.A.S.,  since  pilot  conments  indicated  that  the  authority 
of  the  roll  S.A.S.  was  so  small  that  it  was  essentially  ineffective. 

13.2  Comparison  with  VOLAR  results 

The  VOLAR  simulations  for  the  20  knot  and  35  knot  W.O.D.  conditions  are 
directly  compared  on  Figure  18,  which  shows  the  longitudinal  state  variables 
computed  using  the  optimal  pilot  model.  Figure  20  shows  the  lateral  state 
variables  computed  using  the  classical  pilot  model. 

It  is  evident  from  these  figures  that  the  pilot  employs  larger  random 
control  deflections  for  the  35  knot  condition.  Despite  this,  his  control 
of  the  airplane  is  less  precise.  For  example  from  Figure  18(h)  the  standard 
deviation  of  the  airplane  height  relative  to  the  deck  near  the  instant  of 
touchdown  is  0.6m  at  35  knots  W.O.D.  compared  to  less  than  0.1m  at  20  knots 
W.O.D.  The  pilot's  random  control  activity  has  increased  by  a factor  of 
approximately  6,  as  shown  in  Figure  18(d).  In  addition,  the  simulation  of 
the  lateral  degrees  of  freedom  for  35  knots  W.O.D.  (see  Figure  10)  indicated 
the  lateral  position  error  has  a large  random  component  with  a standard 
deviation  of  2.7  meters,  even  without  including  ship  motion  in  the  simulation. 
The  trends  computed  by  VOLAR  are  thus  in  agreement  with  the  flight  data. 

The  above  correlation  is  encouraging,  but  it  is  clearly  limited  by  the 
fact  that  the  dynamic  interface  chart  is  more  indicative  of  stationkeeping 
performance  than  of  the  complete  recovery  process.  Definitive  criteria  for 
establishing  safe  operational  boundaries  are  not  known.  Detailed  flight 
test  records  are  required  to  discover  the  quantitative  limits  on  state  and 
control  variables  that  determine  the  boundary  between  acceptable  and  unacceptabl 
approaches  and  for  comparison  to  the  simulated  data. 

♦Subsequent  AV-8A's  incorporated  a 3-axis  S.A.S.  similar  to  the  AV-8B 

S.A.S.  described  in  Reference  14,  and  also  employed  a head-up-display. 


14.0  CONCLUSIONS 


The  VOLAR  digital  computer  simulation  program  models  the  dynamics 
of  launch  and  recovery  of  aircraft  operating  from  ships.  The  program 
includes  mathematical  models  for:  the  airframe  and  its  control  system, 
ship  motion,  the  airwake,  atmospheric  turbulence,  the  display  system  and 
the  human  pilot.  Helicopters  or  fixed-wing  aircraft  can  be  simulated. 

The  current  VOLAR  ship  motion  airwake  models  have  been  specially  developed 
for  small  ships  but  the  program  can  simulate  any  type  of  aircraft  and 
ship. 

Using  inputted  statistical  data  on  the  airwake,  ship  motion  and 
other  disturbances,  each  VOLAR  run  computes  time  histories  of  the  means 
and  variances  of  all  system  state  variables.  The  run  time  required  is 
typically  7 percent  of  that  required  for  a 200-run  Monte  Carlo  simulation. 

VOLAR  simulations  of  AV-8A  recoveries  on  a small  LPD-class  ship  show 
good  agreement  with  trends  reported  from  flight  experience.  Variations  of 
Sea  State  parameters  and  wind-over-deck  were  included  in  these  simulations. 

The  VOLAR  simulations  presented  here  demonstrated  the  use  of  two 
alternative  human  pilot  models.  Both  models  are  based  on  published  data. 

One,  the  "classical"  model  involves  a selection  of  parameters  based  on 
published  verbal  rules.  It  is  therefore  dependent  to  some  extent  on  the 
user's  interpretation  of  these  rules.  The  alternative  model,  the  "optimal 
pilot"  model,  is  defined  mathematically.  This  model  is  recommended  because 
it  is  not  subject  to  individual  differences  of  interpretation.  Currently, 
simulations  using  the  optimal  pilot  model  require  approximately  four  times 
the  run  time  required  by  the  classical  model.  Simplification  of  the  optimal 
pilot  model  should  be  explored  in  future  work. 

The  work  reported  here  was  limited  to  developing  the  program  and 
demonstrating  its  use  by  a few  runs  whicn  were  compared  with  AV-8A  flight 
data.  Further  work  should  be  directed  at  using  the  program  to  explore  the 
effects  on  accuracy  of  variations  in  airframe  characteristics,  control  power, 
approach  direction  and  glideslope,  ship  size,  displays,  and  piloting 
techniques. 
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